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CHAPTER  1 

ELEMENTS  OF  GROUND  MOTION 


A.  INTRODUCTION 

Pre-DICE  THROW  II  is  a two-event  series  of  high  explosive  tests 
performed  at  the  White  Sands  Missile  Range  as  a calibration  for  the 
main  DICE  THROW  event.  The  purpose  of  the  tests  was  to  provide  base- 
line data  on  the  site  geology  (cratering  and  ground  motion  character- 
istics) and  to  compare  results  obtained  using  different  explosives 
and  charge  configurations.  The  first  event  consisted  of  a lOO-ton 
surface  tangent  sphere  of  TNT  detonated  by  a single  centrally  located 
booster  charge.  It  was  fired  on  12  August,  1975,  and  will  be  desig- 
nated "Pre-DICE  THROW  II-l".  The  second  consisted  of  a 120-ton  surface 
tangent  domed  cylinder  (length  to  diameter  ratio  0.75)  of  ANFO  (ammo- 
nium nitrate  fuel  oil  mixture).  Designation  of  the  second  event, 
fired  on  22  September  of  the  same  year  and  detonated  by  seven  boosters 
located  at  equal  intervals  along  the  vertical  axis,  will  be  Pre-DICE 
THROW  II-2.  The  operation  was  conducted  at  Queen  15  site  at  the  White 
Sands  Missile  Range,  New  Mexico,  under  the  direction  of  Field  Command, 
Defense  Nuclear  Agency. 

Twenty-five  experimenters  and  support  agencies  participated  in 
the  Pre-DICE  THROW  II  program.  The  Preliminary  Results  Report  (POR 
6904)  was  published  in  September,  1976,  and  describes  the  entire  oper- 
ation, giving  some  test  data  in  preliminary  form.  The  Test  Execution 
Report,  (POR  6965),  in  preparation,  will  describe  the  as-built  pre-test 
experiments  in  detail.  Ground  motion  measurements  will  be  published 
in  a separate  report  by  Waterways  Experiment  Station  (WES),  as  will 
the  results  of  all  experiments  by  other  agencies  and  laboratories. 

The  purpose  of  this  report  is  to  compare  certain  aspects  (such 
as  peak  values)  of  the  ground  motion  data  with  predictions  made  of 
these  aspects  before  the  tests.  Four  agencies  were  charged  with  pro- 
viding predictions:  the  Air  Force  Weapons  Laboratory  (AFWL);  Waterways 
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Experiment  Station  (WES);  R and  0 Associates  (RDA);  and  Field  Command, 
DNA  (FC).  Making  accurate  predictions  is  essential  to  any  test  pro- 
gram in  order  to  provide  correct  gage  ranges  so  that  useful  data  may 
be  obtained,  and  prediction  evaluation  is  a necessary  function. 

Ground  motion  is  a very  complicated  phenomenon  for  which  there  does 
not  exist  good  quantitative  understanding.  Consequently,  prediction 
evaluation  becomes  especially  significant  in  this  area. 

This  report  consists  of  seven  chapters.  The  first  chapter  in- 
cludes a brief  description  of  the  initiation  of  ground  motion  under  a 
high  explosive  detonation.  The  second  chapter  describes  the  ground 
motion  transducers  used  in  Pre-DICE  THROW  II  and  their  outputs.  Chapter 
Three  discusses  the  prediction  methods.  Chapter  Four  presents  a com- 
plete summary  of  data  from  Event  1 and  compares  the  data  with  pre- 
dictions; Chapter  Five  presents  the  same  information  for  Event  II. 

In  Chapter  Six,  data  from  Events  I and  II  are  compared.  Finally, 

Chapter  Seven  gives  a short  resumd  of  results  and  conclusions. 

B.  NATURE  OF  HE  EXPLOSIONS 

Different  types  of  explosives  have  different  physical  detonation 
characteristics.  In  practical  applications,  tlie  configuration  and 
mode  of  initiation  can  also  be  important.  Therefore,  although  the 
two  pre-DICE  THROW  II  charges,  100  tons  of  TNT  and  120  tons  of  ANFO, 
had  about  the  same  energy  yield,  they  may  have  produced  somewhat  dif- 
ferent results  on  and  in  the  ground.  It  is  difficult  to  predict  the 
differences.  Gross  (far  field)  effects  are,  of  course,  yield  dependent 
and  hence  relatively  insensitive  to  charge  configuration.  This  point 
is  discussed  more  fully  in  Chapter  6. 

For  both  Pre-DICE  THROW  II  events,  the  water  table  lay  about  2 
meters  below  the  surface.  The  material  properties  of  the  soil  dif- 
fered slightly  in  the  imo  cases,  but  the  ground  medium  consisced  pri- 
marily of  horizontal  layers  of  silt,  clay,  and  sandy  gravel  down  to 
150  meters,  the  greatest  depth  explored.  Layer  thicknesses  were  of 
the  order  of  tens  of  centimeters. 
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Craters  in  the  two  cases  were  alike  and  differed  from  the  usual 
crater  formation  by  being  broad  and  shallow.  Also,  the  lips  of  the 
craters  appeared  to  be  made  up  of  relatively  intact  layers  folded 
over.  These  craters  give  the  impression  that  near-surface  spalling 
occurred,  which  implies  that  certain  rather  specific  conditions  pertain. 
Development  of  a low  pressure  region  behind  the  air  shock  wave,  coupled 
with  the  bouyant  motion  of  the  fireball,  can  elevate  loose  material  in 
the  spalled  layer,  leading  to  the  well-known  crown  of  jets  and  subse- 
quent throw-out  of  spalled  material. 

The  effects  which  produce  the  crater  demonstrate  by  their  violence 
that  - very  strong  wave  is  propagated  into  the  ground  directly  from 
the  explosive  region.  This  wave  is  often  called  the  "crater-induced" 
ground  wave,  or  more  correctly  the  "direct  ground  wave".  When  the 
explosive  charge  is  in  contact  with  the  ground,  the  ground  wave  is 
initiated  simultaneously  with  the  air  wave. 

Certain  differences  between  the  craters  of  Events  1 and  2 suggest 
that  different  modes  of  coupling  the  ground  wave  into  the  ground 
medium  were  predominant  in  the  two  cases.  In  particular,  the  flat 
bottoni  of  the  Event  1 crater  was  about  2 meters  below  the  water  table, 
While  the  flat  bottom  of  the  Event  2 crater  was  just  ^ the  water  table. 
The  presence  of  the  water  table  and  the  influence  it  had  on  crater  for- 
mation was  significant,  but  it  alone  does  not  determine  the  level  of 
spallation.  Examination  of  the  early  fireball  photographs  (POR  6917) 
shows  that  the  area  of  the  fireball  in  contact  with  the  ground  grows 
relatively  faster  for  Event  2 than  for  Event  1,  suggesting  that  tne 
crater  for  Event  2 may  have  the  larger  diameter,  as  it  does. 

Even  though  our  understanding  of  how  the  crater  is  formed  is 
incomplete,  it  is  obvious  that  a powerful  compression  wave  is  sent 
into  the  ground  at  the  same  time  the  blast  wave  is  propagated  in  air. 

In  air,  which  is  very  compressible,  a shock  wave  forms  which  moves 
supersonically.  In  the  ground,  which  is  relatively  incompressible,  a 
shock  wave  typically  does  not  form  (or  if  it  does,  it  quickly  dies  out). 
The  propagation  of  energy  in  the  ground  is  by  elastic  and  plastic 
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waves,  and  is  much  more  complex  in  nature  than  the  propagation  of  the 
shock  in  air. 

The  movement  of  the  blast  (shock)  wave  across  the  surface  com- 
presses the  ground  and,  if  the  speed  of  the  air  wave  exc  eds  that  of 
an  acoustic  wave  in  the  ground,  then  another  compression  wave  will  be 
produced  in  the  ground.  This  wave  is  called  the  "airblast  induced 
ground  wave".  As  long  as  the  necessary  conditions  exist  at  the  sur- 
face, it  will  be  produced  in  the  ground. 

As  the  airblast  wave  advances,  it  weakens  and  slows  down  (because 
it  is  a shock).  Elastic  ground  waves  travel  at  constant  speed.  There- 
fore, the  ground  wave  eventually  overtakes  and  outruns  the  air  wave. 

The  radius  of  overtaking  is  called  the  "crossover",  and  the  area  beyond 
the  crossover  is  called  the  "outrunning  region". 

There  are  several  other  waves  produced  in  the  groutid:  shear  waves 
and  surface  waves  exist,  as  do  multiple  reflections  of  compression 
waves.  They  all  show  up  on  the  ground  motion  records.  Howevar,  for 
present  purposes,  they  will  be  neglected.  They  do  not  seem  to  play 
significant  roles  in  close-in  ground  motion  phenomenology. 

C.  WAVE  CONFIGURATION  IN  THE  GROUND 

It  is  clear  from  the  preceding  discussion  that  the  ground  motion 
can  become  very  complex.  However,  the  objectives  of  the  Pre-DICE  THROW  II 
ground  motion  experiments  were  to  measure  accurately  the  relationships 
between  the  main  airblast  induced  ground  wave  and  the  main  direct 
ground  wave.  Thus  a very  much  simplified  situation  can  be  depicted. 

The  first  drawing  (Figure  1.1)  shows  an  idealized  situation  in 
which  the  ground  is  a homogeneous  half-space.  Obviously,  as  long  as 
the  airblast  wave  travels  faste*^  than  the  ground  wave,  there  will  be  a 
depth  above  which  an  instrument  will  record  first  the  airblast  induced 
wave  and  then  the  ground  wave.  At  a greatei^  depth,  the  relative  posi- 
tions of  the  two  waves  on  the  record  will  be  reversed.  When  the  direct 
ground  wave  leads  the  airblast  wave,  then  there  cannot  be  an  airblast 
induced  ground  wave  because  the  airblast  wave  is  no  longer  supersonic 
relative  to  the  ground. 
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Figure  1.1.  Configuration  and  Definitions 


The  second  drawing  (Figure  1.2)  shows  the  effect  of  introducing 
a surface  layer  in  which  the  v/ave  propagation  velocity  is  low.  Here 
the  situation  is  more  complex.  In  particular,  it  now  becomes  possible 
for  an  instrument  to  observe  both  air-  and  ground-induced  waves  in  the 
low  velocity  layer;  and  their  relative  positions  can  be  exchanged 
either  by  locating  two  instruments  at  different  depths  and  the  same 
range,  or  at  different  ranges  and  the  same  depth.  Some  examples  will 
be  given  in  the  next  section. 

Note,  too,  that  there  is  a range,  in  the  case  of  the  low  velocity 
surface  layer,  where  the  direct  ground  wave  can  outdistance  the  air 
wave  (Figure  1.21,  but  there  will  still  be  an  airblast  induced  ground 
wave.  Since  low  velocity  surface  layers  are  the  coimion  in  geological 
profiles  of  interest,  the  situation  depicted  in  Figure  1.2  is  a familiar 
one  and  the  cross-over  range  is  a characteristic  of  most  explosions. 

The  range  at  whi-'h  crossover  occurs  can  often  be  calculated  quite 
accurately. 

It  is  clear  from  Figure  1.2  that  the  air-  and  ground- induced  wave 
fronts  always  have  opposite  inclinations  for  the  case  of  a slow  layer 
ever  a faster  layer.  The  radial  (or  horizontal)  components  of  these 
wave  motions  are  generally  directed  outward  and  are  difficult  to  dis- 
tinguish. The  vertical  components,  on  the  other  hard,  are  oppositely 
directed:  the  air-iiicioiced  motion  is  always  downward  and  the  ground- 
induced  motion  is  always  upward.  Therefore,  the  vertical  components 
immediately  identify  the  relative  arrival  times  of  the  two  induced 
waves . 

D.  INSTRUMENT  RESPONSES 

Instruments  (particularly  velocity  gages),  when  placed  at  any  of 
the  positions  indicated  in  Figures  1.1  and  1.2,  will  respond  to  tlie 
vector  sum  of  a’  wave  motions.  Velocity  gages  made  to  respond  only 
to  motion  alon;:  their  principal  axes  will,  of  course,  give  records 
indicating  the  algebraic  addition  of  wave  components  along  those  axes. 
Consequently  the  vertical  component  gage  will  show  strong  positive  and 
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AIR  WAVES 


Figure  1.2.  Effect  of  a Slow  Surface  Layer 


negative  excursions,  while  the  radial  (horizontal)  component  gage  will 
show  mostly  positive  excursions.  Since  the  real  situation  consists 
of  many  wavelets  of  different  types  moving  in  many  directions,  the  gage 
record  can  be  very  ccmolicated. 

In  an  effort  to  simplify  the  picture  but  still  retain  essentials. 
Figure  1.3  shows  schematically  the  algebraic  summation  of  the  vertical 
and  horizontal  components  of  both  the  air-  and  ground-induced  waves  at 
one  of  the  instrument  positions  of  Figure  1.2.  The  heavy  solid  line 
represents,  in  each  case,  the  algebraic  resultant  of  the  two  vertical 
and  two  horizontal  components,  respectively. 

When  reducing  the  data,  the  maxima,  minima,  and  inflection  points 
are  read  off  the  gage  record.  Notice  that  they  do  not  represent  true 
peak  values.  Some  effort  is  made  to  compensate  for  the  unavoidable 
error  by  reading  excursions  from  the  previous  inflection  point  rather 
than  from  the  baseline.  Nevertheless,  there  are  inherent  in  the  data 
reduction  techniques  some  defects  which  tend  to  underestimate  the  ver- 
tical component  peaks;  and  first  to  underestimate,  then  overestimate 
the  radial  component  peaks.  Consequently,  theoretical  calculations 
which  determine  the  induced  wave  forms  separately  may  be  accurate  yet 
not  compare  well  with  the  experiment.  And  by  the  same  argument,  good 
agreement  between  prediction  and  measurement  is  not  necessarily  con- 
firmation of  the  prediction  technique  because  of  problems  of  superposition 

Figure  1.4  illustrates  the  principles  shown  in  Figures  1.2  and  1.3. 

At  the  same  depth  (namely  0.5  meter)  but  at  two  different  ranges,  the 
air-  and  ground-induced  vertical  components  of  velocity  change  places. 

The  interchange  cannot  be  detected  on  the  radial  (horizontal)  wave- 
forms, although  the  two  waves  are  observable.  These  waveforms  are 
traced  from  Figures  4.49,  4.50,  4.61  and  4.62  of  Chapter  4. 

Referring  to  Figure  1.2  once  more,  it  is  obvious  that  at  a given 
range,  the  same  sequence  of  change  in  the  vertical  component  can  be 
found  in  an  array  of  instruments  arranged  vertically.  Figure  1.5 
illustrates  the  vertical  array  concept.  These  waveforms  were  traced 
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from  Figures  4.51,  4.53  and  4.55  of  Chapter  4.  In  the  upper  wave  form 
(a,  at  0.5  meter  depth),  the  air-induced  ground  wave  predominates.  In 
wave  form  b (at  3.7  m depth),  the  direct  induced  ground  wave  is  begin- 
ning to  lead.  In  wave  form  c (at  6.1  m depth),  the  direct  induced 
ground  wave  now  predominates;  and  in  wave  form  d (at  9.1  m depth),  the 
air-induced  ground  wave  has  essentially  disappeared.  Perhaps  the 
boundary  between  the  slow  and  fast  media  lies  at  or  above  the  9.1 
meter  level. 

The  features  which  have  just  been  described  are  mainly  character- 
istics of  the  slow  layer.  As  long  as  the  induced  waves  form  an  "X"  in 
the  slow  layer  as  shown  in  Figure  1.2,  these  features  will  be  seen. 
Eventually  the  air  wave  falls  behind  the  ground  wave  and  the  "X"  can 
no  longer  form. 

The  air-induced  wave  in  the  slow  layer  must,  of  course,  intersect 
the  slow-fast  interface.  This  was  omitted  from  Figure  1.2  for  sim- 
plicity. At  the  slow-fast  interface,  the  air-induced  wave  will  be  both 
reflected  and  refracted.  The  refracted  wave  will  penetrate  the  fast 
layer,  but  the  laws  of  elastic  refraction  dictate  that  the  transmitted 
wave  be  relatively  weak.  Therefore,  the  slow  layer  tends  to  isolate 
the  region  beneath  it  from  the  air- induced  wave,  and  this  can  be  seen 
by  examination  of  Figure  1.5. 

E.  RELATIONS  TO  PREDICTIONS 

Four  prediction  techniques  will  be  discussed  in  some  detail  in 
Chapter  3.  Here  we  deal  in  generalities. 

Predictions  which  separately  calculate  the  air-  and  direct-induced 
wave  shapes  by  hydrocode  or  some  other  method  based  on  first  principles, 
but  do  not  combine  them  in  proper  phase,  can  be  expected  to  overpredict 
peak  values  for  vertical  velocity  components  and  to  underpredict  peak 
values  for  radial  (horizontal)  velocity  components.  This  is  illustrated 
in  Figure  1 .3. 

If  such  predictions  include  combining  the  air-  and  direct-induced 
wave  shapes  in  proper  phase,  then  the  predictions  will  be  very  sensitive 
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to  the  integral  of  the  propagation  velocities  over  the  full  ray  paths 
of  both  the  airblast  wave  and  the  direct  ground  wave. 

Prediction  techniques  based  on  results  from  other  tests  may  give 
excellent  or  absurd  results  depending  on  charge  configuration  and  geo- 
logic similarity  considerations. 

It  will  be  shown  in  Chapters  4 and  5 that  predictions  of  phase 
relations  between  air-  and  direct-induced  waves  were  not  very  success- 
ful. Predictions  of  peak  values  versus  distance  were  in  fairly  good 
agreement  with  measured  peak  values  when  phase  was  ignored.  The  first 
statement  attests  to  the  difficulty  of  integrating  along  the  ray  path. 
The  second  statement  attests  to  the  generally  similar  nature  of  these 
field  tests  to  others  performed  over  media  of  approximately  the  same 
geologic  profile. 
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CHAPTER  2 


GROUND  MOTION  DATA  ACQUISITION  TECHNIQUE 

A.  INTRODUCTION 

The  first  event  (TNT)  on  Pre-DICE  THROW  II  had  the  most  ground 
motion  instrumentation.  There  were  a total  of  78  channels  of  accelera- 
tion measurements  and  110  channels  of  velocity  measurements.  The  sec- 
ond event  (ANFO)  had  32  channels  of  acceleration  and  72  channels  of 
velocity  measurements.  A break-out  of  gage  placement  for  each  event 
is  shown  in  Tables  2.1  through  2.4.  Figures  2.1  and  2.2  indicate  the 
placement  pattern  surrounding  ground  zero  for  each  event  and  Figures 
2.3  and  2.4  depict  the  gage  placement  profiles  (range  and  depth)  for 
each  event. 

Additional  radial s (60,  180,  and  300  degrees)  were  instrumented 
on  the  TNT  event  at  the  33.5-meter  range  to  examine  test-bed  azimuthal 
symmetry. 

B.  TRANSDUCERS  AND  THEIR  SYSTEMS 

1 . Accelerometers 

Endevco  piezoresisti ve  accelerometers  were  used  on  both 
events.  The  different  models  used  were  chosen  for  the  following  reasons: 

a.  The  gage's  range  is  selected  so  that  it  is  4 to  10  times 
greater  than  the  set  range  or  predicted  level  of  acceleration.  This  is 
done  for  several  reasons:  (1)  the  h'  ' >r  the  gage's  range,  the  higher 
is  its  frequency  response;  (2)  above  ^»ercent  of  the  gage's  range, 
the  outputs  are  sufficient  to  give  a good  signal  to  noise  ratio;  and 
(3)  additional  insurance  is  provided  against  overranging. 

b.  Due  to  lead-time  requirements  by  the  manufacturer,  the 
order  for  the  gages  was  placed  based  on  preliminary  ground  motion  pre- 
dictions. Up-dating  of  these  predictions  caused,  in  some  cases,  re- 
arrangement of  the  gages  to  different  locations  more  appropriate  to 
their  range. 
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c.  The  miniature  model  2264  gages  were  used  in  the  close-in 
region  because  of  their  high  natural  frequency.  These  gages  were  houseu 
in  micro  canisters,  which  increases  the  overall  natural  frequency  of 
the  system. 

The  overall  system  schematic  used  for  accelerometer  channels  is 
shown  in  Figure  2.5.  The  data  analysis  techniques  will  be  discussed 
later. 

2.  Velocity  Gages 

The  velocity  gages  used  were  versions  of  the  Sandia  DX  gage 
and  were  manufactured  by  Consolidated  Electrodynamics  Corporation  and 
Sparton  Southwest,  Inc.  Only  two  types  of  these  gages  were  used.  In 
general,  the  close-in  gages  had  smaller  sensing  armatures,  thereby 
giving  them  a lower  sensitivity.  All  of  the  gages  were  damped  using 
1000  cs  damping  fluid.  The  system  schematic,  shown  in  Figure  2.5, 
also  applies  to  the  velocity  measurements. 

3.  Canisters 

The  three  types  of  canisters  used  by  WES  for  housing  the 
ground  motion  transducers  are  called  "micro",  "mini"  and  "standard". 

The  "micro"  and  "mini"  canisters  are  used  exclusively  for  housing 
accelerometers,  and  the  "standard"  canister  can  contain  two  velocity 
gages  and  two  accelerometers.  The  fill  or  potting  material  used  inside 
the  canisters  is  a polyurethane-based  thermosetting  resin  called  Siwax 
601.  A design  and  evaluation  discussion  on  these  canisters  is  found  in 
Reference  1 . 

4.  Grout 

All  ground  motion  canisters  were  fixed  in  place  in  boreholes 
with  a "quick-set"  grout.  A plug  of  this  grout,  roughly  0.61  m long, 
and  extending  from  about  0.30  m below  canister  center  to  0.30  m above 
it,  was  used.  A filler  grout,  designated  WES  grout  E-2-E,  was  used 


^Ground  Motion  Canister  Design  and  Evaluation,  Andres  Peekna,  U.S.  Army 
Engineer  Waterways  Experiment  Station  Corps  of  Engineers,  Vicksburg, 
Mississippi . 
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between  gages,  and  to  top  out  deep  (3.7  m and  greater)  canister  holes. 
Shallow  (1.8  m and  less)  holes  were  topped  out  with  local  top-soil 
backfill.  Some  properties  of  the  grouts  mentioned  are: 


Density 


Quick  Set 
2026.3  kg/m3 


E-2-E 

2023.1  kg/m3 


Unconfined  13,800-20,700  kPa  690-1035  kPa 

Compressive 

Strength 


Mix  Components 
(for  0.30  m^ 
batch) 


Chem  Stress  Cement  5.7  kg 
Cal -Seal  Cement  17.2  kg 
Sand  22.9  kg 
Water  11.5  kg 


Portland  Cement  5.6  kg 
Bentolite  Gel  2.1  kg 
Barite  8.3  kg 
Sand  28.5  kg 
Water  12.8  kg 


5 . Gage  Emplacement 

A special  placement  rod  was  used  to  lower  the  gage  canisters 
into  the  holes.  Each  section  of  this  rod  was  0.30  m in  length.  The 
rods  and  canister  had  alignment  marks  which  were  aligned  throughout 
the  lowering  process.  When  the  canister  was  at  the  proper  depth,  a 
4-power  scope  was  attached  to  the  top  rod  and  alignment  with  ground 
zero  was  made.  This  alignment  was  maintained  during  the  grouting  process. 

The  accelerometer  canisters  were  aligned  vertically  and  hori- 
zontally using  the  assumption  that  the  holes  were  vertical  (WES  feels 
that  their  canisters  were  within  2 to  3 degrees  from  vertical).  The 
errors  in  gage  output  due  to  alignment  are  negligible  if  this  degree 
of  alignment  accuracy  is  maintained. 

Since  the  velocity  gage  is  extremely  sensitive  with  respect 
to  horizontal  and  vertical  orientations,  a bubble  table  was  used  prior 
to  canister  emplacement.  The  canister  was  mounted  on  the  table,  and 
the  electrical  zero  of  each  gage  was  determined.  After  the  canisters 
were  lowered,  this  same  electrical  zero  was  maintained  during  grouting. 

The  deviations  from  true  vertical  or  horizontal  were  maintained  at  less 
than  1/4  degree. 


2-3 


C.  METHOD  USED  BY  WES  TO  REDUCE  DATA 


1 . Accelerometer  Data 

On  horizontal  acceleration  data,  the  outward  peak  (away  from 
ground  zero)  was  normally  taken  as  the  peak  acceleration.  For  ver- 
tical acceleration,  the  absolute  value  of  the  largest  peak  in  either 
direction  was  taken  as  the  peak  value. 

The  raw  data  was  digitized  at  the  rate  of  24,000  samples  per 
second.  No  filtering  was  used. 

2.  Velocity  Gage  Data 

The  outward  peak  motion  was  chosen  for  the  horizontal  velocity 
peak  value.  The  inward  peak  (toward  ground  zero)  was  considered  as  an 
elastic  recovery  of  material  and  was  normally  less  than  the  outward 
excursion.  Two  values  were  read  from  the  vertical  velocity  records, 
the  downward  directed  airblast  induced  peak  and  the  positive  peak  caused 
by  direct  or  crater  induced  ground  motion.  Selecting  or  interpreting 
these  values  prior  to  the  outrunning  region  was  relatively  straight- 
forward; however,  for  those  data  at  ranges  just  past  outrunning  the 
task  was  more  difficult,  sometimes  requiring  large  expansions  (in  time) 
of  the  waveform  in  order  to  discern  the  change  in  slope  cf  the  data. 

The  values  chosen  for  peaks  are  not  necessarily  measured  from  zero,  but 
are  taken  from  the  point  where  the  airblast  induced  velocity  starts 
to  influence  the  waveshape  caused  by  the  ground  motion  induced  portion. 
As  discussed  in  Chapter  1,  some  error  in  inevitable  using  i..iis  method 
if  direct  and  airblast  induced  motions  are  superimposed. 

3.  Record  Quality 

Peak  data  was  not  obtained  from  21  accelerometers  on  Event  1, 
and  5 accelerometers  on  Event  2.  In  most  cases  this  was  because  the 
cables  or  cable/transducer  interfaces  were  destroyed  prior  to  peak 
arrival.  Band  edge  was  exceeded  on  16  records  on  Event  1 and  12  rec- 
ords on  Event  2.  Table  2.5  gives  specific  locations  of  those  gages 
from  which  band  edge  data  or  records  indicating  cable  destruction  were 
obtained.  Note  that  all  gages  in  the  band  edge  category  are  below  the 
water  table. 
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Six  vertical  velocity  gages  went  to  band  edge  on  Pre-DICE 
THROW  II,  Event  1.  Five  of  these  were  at  the  48.8-m  range  on  all  three 
main  gage  radials.  They  went  out  of  band  edge  on  the  negative  peaks 
(airblast  induced)  only  and  subsequent  portions  of  the  wave  forms  were 
within  band  edge.  The  depth  of  these  gages  was  0.5  to  1.8  meters. 

The  signal  conditioning  for  similar  gages  on  Event  2 was 
changed  so  that  the  gage  overranging  would  not  occur. 


Table  2.1.  Pre-DICE  THROW  II,  Event  1 
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ACCELEROMETER  IDENTIFICATION 


Radial 

Range 

Depth 

Orientation 

Gage 

Gage 

Set 

Range 

Canister 

(deg) 

(m) 

(m) 

Mode* 

Model** 

Range 

(g's) 

Type 

0 

12.2 

0.5 

V 

2264 

30  kg 

10  kg 

Micro 

L. 

0 

12.2 

0.5 

H 

2264 

20  kg 

5.5  kg 

Micro 

P 

i 0 

12.2 

1.8 

V 

2264 

20  kg 

10  kg 

Micro 

f 

0 

12.2 

1.8 

H 

2264 

20  kg 

4.5  kg 

Micro 

e 

i 

0 

12.2 

3.7 

V 

2264 

30  kg 

10  kg 

Micro 

1 

1 

0 

; 

12.2 

3.7 

H 

2264 

20  kg 

4.5  kg 

Micro 

13 

0 

18.3 

0.5 

V 

2264 

5 kg 

2.3  kg 

Mi  cro 

1 

\ 

0 

18.3 

0.5 

H 

2264 

5 kg 

1.3  kg 

Micro 

f 

\ 

0 

18.3 

1.8 

V 

2264 

5 kg 

2 kg 

Micro 

0 

18.3 

1.8 

H 

2264 

5 kg 

1.1  kg 

Micro 

1 

0 

18.3 

3.7 

V 

2261M6 

10  kg 

2 kg 

Mini 

''r 

0 

18.3 

3.7 

H 

2261M6 

10  kg 

1.1  kg 

Mini 

: 

0 

18.3 

6.1 

V 

2261M6 

10  kg 

2 kg 

Mini 

' 

0 

18.3 

6.1 

H 

2261M6 

10  kg 

1.1  kg 

Mini 

f- 

L 

0 

18.3 

9.1 

V 

2261M6 

10  kg 

2 kg 

Mini 

0 

18.3 

9.1 

H 

2261M6 

10  kg 

1.1  kg 

Mini 

0 

9.1 

3.7 

V 

2261C 

2.5  kg 

700  g 

Mini 

0 

9.1 

3.7 

H 

2261 C 

2.5  kg 

400  g 

Mini 

L 

0 

24.4 

6.1 

V 

2261C 

2.5  kg 

700  g 

Mini 

0 

24.4 

6.1 

H 

2261 C 

2,5  kg 

400  g 

Mini 

} 

0 

24.4 

9.1 

V 

2261C 

2.5  kg 

700  g 

Mini 

1 

0 

24.4 

9.1 

H 

2261C 

2.5  kg 

400  g 

Mini 

1 

0 

24.4 

12.2 

V 

2261C 

2.5  kg 

700  g 

Mini 

0 

24.4 

12.2 

H 

2261C 

2.5  kg 

400  g 

Mini 

‘ 

0 

24.4 

18.3 

V 

2261C 

2.5  kg 

700  g 

Mini 

£ 

0 

24.4 

18.3 

H 

2261C 

2.5  kg 

400  g 

Mini 

i.; 

0 

33.6 

0.5 

V 

2261C 

2.5  kg 

300  g 

Std. 

|: 

0 

33.5 

0.5 

H 

2261C 

2.5  kg 

170  g 

Std. 
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Table  2.1.  Pre-DICE  THROW  II,  Event  1 (Continued) 
ACCELEROMETER  IDENTIFICATION 


Radial 

Range 

Depth 

Orientation 

Gage 

Gage 

Set 

Range 

Canister 

(deg) 

(mj 

(m) 

Mode* 

Model** 

Range 

(g's) 

Type 

0 

33.5 

3.7 

V 

2262L 

1 kg 

200  g 

Std. 

0 

33.5 

3.7 

H 

2262C 

200  g 

150  g 

Std. 

0 

33.5 

6.1 

V 

2261C 

2.5  kg 

200  g 

Mini 

0 

33.5 

6.1 

H 

2262C 

200  g 

150  g 

Mini 

0 

33.5 

9.1 

V 

2261 C 

2.5  kg 

200  g 

Mini 

0 

33.5 

9.1 

H 

2262C 

' 200  g 

150  g 

Mini 

0 

33.5 

12.2 

V 

2261 C 

2.5  kg 

200  g 

Mini 

0 

33.5 

12.2 

H 

2262C 

200  g 

150  g 

Mini 

0 

48.8 

0.5 

V 

2262C 

1 kg 

140  g 

Std. 

0 

48.8 

0.5 

H 

2262C 

200  g 

45  g 

Std. 

0 

48.8 

3.7 

V 

2262C 

200  g 

75  g 

Std. 

i ° 

48.8 

3.7 

H 

2262C 

200  g 

40  g 

Std. 

0 

91.4 

0.5 

V 

2262C 

200  g 

60  g 

Std. 

‘ 0 

91.4 

0.5 

H 

2262C 

200  g 

15  g 

Std. 

r 

0 

91.4 

1.8 

V 

2262C 

200  g 

25  g 

Std. 

0 

91.4 

1.8 

H 

2262C 

25  g 

8 g 

Std. 

0 

91.4 

3.7 

V 

2262C 

200  g 

25  g 

Std. 

0 

91.4 

3.7 

H 

2262C 

25  g 

8 g 

Std. 

0 

91.4 

6.1 

V 

2262C 

200  g 

25  g 

Std. 

0 

91.4 

6.1 

H 

2262C 

25  g 

8 g 

Std. 

60 

33.5 

3.7 

V 

2262C 

1 kg 

200  g 

Mini 

i 60 

33.5 

3.7 

H 

2262C 

200  g 

150  g 

Mini 

; 120 

18.3 

0.5 

V 

2264 

5 kg 

2.3  kg 

Micro 

120 

18.3 

0.5 

H 

2264 

5 kg 

1.3  kg 

Micro 

120 

24.4 

3.7 

V 

2261 C 

2.5  kg 

700  g 

Mini 

120 

24.4 

3.7 

H 

2261C 

2.5  kg 

400  g 

Mini 
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Table  2.1.  Pre-DICE  THROW  II » Event  1 (Continued) 
ACCELERATION  IDENTIFICATION 


Set 


Radial 

(deg) 

Range 

(m) 

Depth 

(m) 

Orientation 

Mode* 

Gage 

Model** 

Gage 

Range 

Range 

(g's) 

Canister 

Type 

120 

24.4 

6.1 

V 

2261 C 

2.5  kg 

700  g 

Mini 

120 

24.4 

6.1 

H 

2261 C 

2.5  kg 

400  g 

Mini 

120 

24.4 

9.1 

V 

2261 C 

2.5  kg 

700  g 

Mini 

120 

24.4 

9.1 

H 

2261 C 

2.5  kg 

400  g 

Mini 

120 

33.5 

3.7 

V 

2262C 

1 kg 

200  g 

Mini 

120 

33.5 

3.7 

H 

2262C 

200  g 

150  g 

Mini 

120 

33.5 

6.1 

V 

2261 C 

2.5  kg 

200  g 

Mini 

120 

33.5 

6.1 

H 

2262C 

200  g 

150  g 

Mini 

180 

33.5 

3.7 

V 

2262C 

1 kg 

200  g 

Mini 

180 

33.5 

3.7 

H 

2262C 

200  g 

150  g 

Miri 

240 

18.3 

0.5 

V 

2264 

5 kg 

2.3  kg 

Micro 

240 

18.3 

0.5 

H 

2264 

5 kg 

1.3  kg 

Micro 

240 

24.4 

3.7 

V 

2261 C 

2.5  kg 

700  g 

Mini 

240 

24.4 

3.7 

H 

2261 C 

2.5  kg 

400  g 

Mini 

240 

24.4 

6.1 

V 

2261C 

2.5  kg 

700  g 

Mini 

240 

24.4 

6.1 

H 

2261 C 

2.5  kg 

400  g 

Mini 

240 

24.4 

9.1 

V 

2261C 

2.5  kg 

700  g 

Mini 

240 

24.4 

9.1 

H 

226 1C 

2.5  kg 

400  g 

Mini 

240 

33.5 

3.7 

V 

2262C 

1 kg 

200  g 

Mini 

240 

33.5 

3.7 

H 

2262C 

200  g 

150  g 

Mini 

240 

33.5 

6.1 

V 

2261 C 

2.5  kg 

200  g 

Mini 

240 

33.5 

6.1 

H 

2262C 

200  g 

150  g 

Mini 

300 

33.5 

3.7 

V 

2202C 

1 kg 

200  g 

Mini 

300 

33.5 

3.7 

H 

2262C 

200  g 

150  g 

Mini 

* H - Horizontal 
V - Vertical 

**  Manufactured  by  Endevco 
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Table  2.2.  Pre-OICE  THROW  II,  Event  1 


velocity  gage  identification 


Radial 

(deg) 

Range 

(m) 

Depth 

(m) 

Orientati 

Mode* 

0 

24.4 

0.5 

V 

0 

24.4 

0.5 

H 

0 

24.4 

0.5 

T 

0 

24.4 

1.8 

V 

0 

24.4 

1.8 

H 

0 

33.5 

0.5 

V 

0 

33.5 

0.5 

H 

0 

33,5 

1.8 

V 

0 

33.5 

1.8 

H 

0 

33.5 

3.7 

V 

0 

33.5 

3.7 

H 

0 

48.8 

0.5 

V 

0 

48.8 

0.5 

H 

0 

48.8 

1.8 

V 

0 

48.8 

1.8 

H 

0 

48.8 

3.7 

V 

0 

48.8 

3.7 

H 

0 

48.8 

6.1 

V 

0 

48.8 

b.l 

H 

0 

48.8 

9.1 

V 

0 

48.8 

9.1 

H 

0 

70.1 

0.5 

V 

0 

70.1 

0.5 

H 

0 

70.1 

0.5 

T 

0 

/O.l 

1.8 

V 

0 

70.1 

1.8 

H 

0 

70.1 

3.7 

V 

Fluid 

Viscosity 

(cs)** 

Gage 

Set 

Range 

(mps) 

Canister 

Type 

L-1 

5.5 

Std. 

L-1 

3.5 

Std. 

L-1 

3.5 

Std. 

L-1 

3.7 

Std. 

L-1 

3.5 

Std. 

L-1 

3.0 

Std. 

L-1 

1.8 

Std. 

L-1 

2.1 

Std. 

L-1 

1.8 

Std. 

L-1 

1.5 

Std. 

L-1 

1.8 

Std. 

H-1 

1.5 

Std. 

H-1 

0.82 

Std. 

H-1 

1.1 

Std. 

H-1 

0.82 

Std. 

H-1 

0.76 

Std. 

H-1 

0.82 

Std. 

H-1 

0.49 

Std. 

H-1 

0.82 

Std. 

H-1 

0.30 

Std. 

H-1 

0.82 

Std. 

H-1 

0.76 

Std. 

H-1 

0.43 

Std. 

H-1 

0.43 

Std. 

H-1 

0.55 

Std. 

H-1 

0.43 

Std. 

H-1 

0.40 

Std. 

2-9 


Table  2.2.  Pre-DICE  THROW  II,  Event  1 (Continued) 
VELOCITY  GAGE  IDENTIFICATION 


Radial 

(deg) 

Range 

(m) 

Deptn 

(m) 

Orientation 

Mode* 

Fluid 

Viscosity 

(cs)** 

Gage 

Set 

Range 

(mps) 

Canister 

Type 

0 

70.1 

3.7 

V 

H-1 

0.40 

Std. 

0 

70.1 

3.7 

H 

H-1 

0.43 

Std. 

0 

70.1 

6.1 

V 

H-1 

0.27 

Std. 

0 

70.1 

6.1 

H 

H-1 

0.43 

Std. 

0 

91.4 

0.5 

V 

H-1 

0.49 

Std. 

0 

91.4 

0.5 

H 

H-1 

0.21 

Std. 

0 

91.4 

1.8 

V 

H-1 

0.37 

Std. 

0 

91.4 

1.8 

H 

H-1 

0.21 

Std. 

0 

91.4 

3.7 

V 

H-1 

0.26 

Std. 

0 

91.4 

3.7 

H 

H-1 

0.21 

Std. 

0 

91.4 

6.1 

V 

H-1 

0.17 

Std. 

0 

91.4 

6.1 

H 

H-1 

0.21 

Std. 

0 

121.9 

0.5 

V 

H-1 

0.40 

Std. 

0 

121.9 

0.5 

H 

H-1 

0.18 

Std. 

0 

121.9 

1.8 

V 

H-1 

0.30 

Std. 

0 

121.9 

1.8 

H 

H-1 

0.14 

Std. 

0 

121.9 

3.7 

V 

H-1 

0.22 

Std. 

0 

121.9 

3.7 

H 

H-1 

0.12 

Std. 

0 

121.9 

6.1 

V 

H-1 

0.15 

Std. 

0 

121.9 

6.1 

H 

H-1 

0.12 

Std. 

0 

121.9 

9.1 

V 

H-1 

0.12 

Std. 

0 

121.9 

9.1 

H 

H-1 

0.12 

Std. 

0 

121.9 

12.2 

V 

H-1 

0.12 

Std. 

0 

121.9 

12.2 

H 

H-1 

0.12 

Std. 

0 

182.9 

0.5 

V 

H-1 

0.34 

Std. 

0 

182.9 

0.5 

H 

H-1 

0.13 

Std. 

0 

182.9 

1.8 

V 

H-1 

0.24 

Std. 
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Table  2.2.  Pre-DICE  THROW  II,  Event  1 (Continued) 

VELOCITY  GAGE  IDENTIFICATION 

Gage 

Fluid  Set 

Range  Depth  Orientation  Viscosity  Range  Canister 


Table  2.2.  Pre-DICE  THROW  II,  Event  1 (Continued) 


VELOCITY  (iAGE  IDENTIFICATION 


Radial 

(deg) 

Range 

(tn) 

Depth 

(m) 

Orientation 

Mode* 

Fluid 

Viscosity 

(cs)** 

Gage 

Set 

Range 

(mps) 

Canister 

Type 

0 

182.9 

1.8 

H 

H-1 

0.10 

Std. 

0 

182.9 

3.7 

V 

H-1 

0.18 

Std. 

0 

182.9 

3.7 

H 

H-1 

0.07 

Std. 

60 

33.5 

0.5 

V 

L-1 

3.0 

Std. 

60 

33.5 

0.5 

H 

L-1 

1.8 

Std. 

120 

24.4 

0.5 

V 

L-1 

5.5 

Std. 

120 

24.4 

0.5 

H 

L-1 

3.5 

Std. 

120 

24.4 

1.8 

V 

L-1 

3.7 

Std. 

120 

24.4 

1.8 

H 

L-1 

3.5 

Std. 

120 

33.5 

0.5 

V 

L-1 

3.0 

Std. 

120 

33.5 

0.5 

H 

L-1 

1.8 

Std. 

120 

33.5 

1.8 

V 

L-1 

2.1 

Std. 

120 

33.5 

1.8 

H 

L-1 

1.8 

Std. 

120 

48.8 

0.5 

V 

H-1 

1.5 

Std. 

120 

48.8 

0.5 

H 

H-1 

0.82 

Std. 

120 

48.8 

1.8 

V 

H-1 

1.1 

Std. 

120 

48.8 

1.8 

H 

H-1 

0.82 

Std. 

120 

48.8 

3.7 

V 

H-1 

0.76 

Std. 

120 

48.8 

3.7 

H 

H-1 

0.82 

Std. 

120 

70.1 

0.5 

V 

H-1 

0.76 

Std. 

120 

70.1 

0.5 

H 

H-1 

0.43 

Std. 

120 

70.1 

3.7 

V 

H-1 

0.40 

Std. 

120 

70.1 

3.7 

H 

H-1 

0.43 

Std. 

120 

121.9 

0.5 

V 

H-1 

0.40 

Std. 

12L- 

121.9 

0.5 

■| 

H-1 

0.18 

Std. 

120 

121.^ 

1.8 

V 

H-1 

0.30 

Std. 

120 

121.9 

1.8 

H 

H-1 

0.14 

Std. 
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Table  2.2.  Pre-DICE  THROW  II,  Event  1 (Continued) 
VELOCITY  GAGE  IDENTIFICATION 


Radial 

(deg) 

Range 

(m) 

Depth 

(m) 

Orientation 

Mode* 

F^uid 

Viscosity 

(cs)** 

Gage 

Set 

Range 

(mps) 

Canister 

Type 

240 

121.9 

6.1 

V 

H-1 

0.15 

Std. 

240 

121.9 

6.1 

H 

H-1 

0.12 

Std. 

300 

33.5 

0.5 

V 

L-1 

3.0 

Std. 

300 

33.5 

0.5 

H 

L-1 

1.8 

Std. 

* H - Horizontal 

V - Vertical 

T - Transverse 

**  H - High  Sensitivity 

L - Low  Sensitivity 

-1  - 1000  cs  damping  fluid 
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Table  2.3.  Pre-DICE  THROW  II,  Event  2 
ACCELEROMETER  IDENTIFICATION 


i 


Radial 

(deg) 

Range 

(m) 

Depth 

(m) 

Orientation 

Mode* 

Gage 

Model** 

Gage 

Range 

(g) 

Set 

Range 

(g) 

Canister 

Type 

0 

12.2 

0.5 

V 

2264 

30  kg 

10  kg 

Micro 

0 

12.2 

0.5 

H 

2264 

20  kg 

5.5  kg 

Micro 

0 

12.2 

1.8 

V 

2264 

30  kg 

10  kg 

Micro 

0 

12.2 

1.8 

H 

2264 

20  kg 

5.5  kg 

Micro 

0 

18.3 

0.5 

V 

2264 

5 kg 

2.3  kg 

Micro 

0 

18.3 

0.5 

H 

2264 

5 kg 

1 .3  kg 

Micro 

0 

18.3 

1.8 

V 

2264 

5 kg 

2 kg 

Micro 

0 

18.3 

1.8 

H 

2264 

5 kg 

1 .1  kg 

Mi  cro 

0 

18.3 

3.7 

V 

2261 M6 

10  kg 

2 kg 

Mini 

0 

18.3 

3.7 

H 

2261M6 

10  kg 

1.1  kg 

Mini 

0 

18.3 

6.1 

V 

2261 M6 

10  kg 

2 kg 

Mini 

0 

18.3 

6.1 

H 

2261M6 

10  kg 

1 .1  kg 

Mini 

0 

18.3 

9.1 

V 

2261M6 

10  kg 

2 kg 

Mini 

0 

18.3 

9.1 

H 

2261M6 

10  kg 

1 .1  kg 

Mini 

0 

24.4 

3.7 

V 

226 1C 

2.5  kg 

700  a 

Mini 

0 

24.4 

3.7 

H 

226 1C 

2.5  kg 

400  g 

Mini 

0 

24.4 

6.1 

V 

2261C 

2.6  kg 

700  g 

Mini 

0 

24.4 

6.1 

H 

2261 C 

2.5  kg 

400  g 

Mini 

0 

24.4 

9.1 

V 

226 1C 

2.5  kg 

700  g 

Mini 

0 

24.4 

9.1 

H 

226 1C 

2.5  kg 

400  g 

Mini 

0 

33.5 

3.7 

V 

2262C 

1 kg 

200  g 

Mini 

0 

33.5 

3.7 

H 

2262C 

200  g 

150  g 

Mini 

0 

33.5 

6.1 

V 

2262C 

1 kg 

200  g 

Mini 

0 

33.5 

6.1 

H 

2262C 

200  g 

150  g 

Mini 

120 

33.5 

3.7 

V 

2262C 

1 kg 

200  g 

Mini 

120 

33.5 

3.7 

H 

2262C 

200  g 

150  g 

Mini 

120 

33.5 

6.1 

V 

2262C 

1 kg 

200  g 

Mini 

f 

-5 
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Table  2.3.  Pre-DICt  THROW  II,  Event  2 (Continued) 
ACCELEROMETER  IDENTIFICATION 


Radial 

(deg) 

Range 

(m) 

Depth 

(m) 

Orientation 

Mode* 

120 

33.5 

6.1 

H 

240 

33.5 

3.7 

V 

240 

33.5 

3.7 

H 

240 

33.5 

6.1 

V 

240 

33.5 

6.1 

H 

* H - Horizontal 
V - Vertical 

**  Manufactured  by  Endevco 


Gage 

Model** 

Gage 

Range 

(g) 

Set 

Range 

(g) 

Canister 

Type 

2262C 

200  g 

150  g 

Mini 

2262C 

1 kg 

200  g 

Mini 

2262C 

200  g 

150  g 

Mini 

2262C 

1 kg 

200  g 

Mini 

:;262C 

200  g 

150  g 

Mini 

> 
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Table  2.4.  Pre-DICE  THROW  II,  Event  2 
VELOCITY  GAGE  IDENTIFICATION 


I 

I 

I 

I 

I 


ij 

: 


Gage 


Radial 

Range 

Depth 

Orientation 

Fluid 

Viscosity 

Set 

Range 

Canister 

(deg) 

(m) 

(m) 

Mode* 

(cs)** 

(mps) 

Type 

0 

24.4 

0.5 

V 

L-1 

5.5 

Std. 

0 

24.4 

0.5 

H 

L-1 

3.5 

Std. 

0 

24.4 

1.8 

V 

i.-l 

3.7 

Std. 

0 

24.4 

1.8 

H 

L-1 

3.5 

Std. 

0 

33.5 

0.5 

V 

L-1 

3.0 

Std. 

0 

33.5 

0.5 

H 

L-1 

1.8 

Std. 

0 

33.5 

1.8 

V 

L-1 

2.1 

Std. 

0 

33.5 

1.8 

H 

L-1 

1.8 

Std. 

0 

33.5 

9.1 

V 

H-1 

0.55 

Std. 

0 

33.5 

9.1 

H 

H-1 

1.1 

Std. 

0 

39.6 

0.5 

V 

L-1 

4.3 

Std. 

0 

39.6 

0.5 

H 

L-1 

2.1 

Std. 

0 

48.8 

0.5 

V 

H-1 

1.5 

Std. 

0 

48.8 

0.5 

H 

H-1 

0.82 

Std. 

0 

48.8 

1.8 

V 

H-1 

1.1 

Std. 

0 

48.8 

1.8 

H 

H-1 

0.82 

Std. 

0 

48.8 

3.7 

V 

H-1 

0.76 

Std. 

0 

48.8 

3.7 

H 

H-1 

0.82 

Std. 

0 

48.8 

6.1 

V 

H-1 

0.49 

Std. 

0 

48.8 

6.1 

H 

H-1 

0.82 

Std. 

0 

48.8 

9.1 

V 

H-1 

0.30 

Std. 

0 

48.8 

9.1 

H 

H-1 

0.82 

Std. 

0 

57.9 

0.5 

V 

L-1 

3.4 

Std. 

0 

57.9 

0.5 

H 

H-1 

0.91 

Std. 

0 

70.1  , 

0.5 

V 

H-1 

0.76 

Std. 

0 

70.1 

0.5 

H 

H-1 

0.43 

Std. 

0 

70.1 

1.8 

V 

H-1 

0.55 

Std. 

0 

70.1 

1.8 

H 

H-1 

0.43 

Std. 

Table  2.4.  Pre-DICE  THROW  II,  Event  2 (Continued) 
VELOCITY  GAGE  IDENTIFICATION 


Radial 

(deg) 

Range 

(>n) 

Depth 

(m) 

Orientation 

Mode* 

Fluid 

Viscosity 

(cs)** 

Gage 

Set 

Range 

(mps) 

Canister 

Type 

0 

70.1 

3.7 

V 

H-1 

0.40 

Std. 

0 

70.1 

3.7 

H 

H-1 

0.43 

Std. 

0 

70.1 

6.1 

V 

H-1 

0.27 

Std. 

0 

70.1 

6.1 

H 

H-1 

0.43 

Std. 

0 

85.3 

0.5 

V 

H-1 

1.1 

Std. 

0 

85.3 

0.5 

H 

H-1 

0.30 

Std. 

0 

91.4 

0.5 

V 

H-1 

0.91 

Std. 

0 

91.4 

0.5 

H 

H-1 

0.21 

Std. 

0 

97.5 

0.5 

V 

H-1 

0.67 

Std. 

0 

97.5 

0.5 

H 

H-1 

1.8 

Std. 

0 

106.7 

0.5 

V 

H-1 

0.55 

Std. 

0 

106.7 

0,5 

H 

H-1 

0.15 

Std. 

0 

121.9 

0.5 

V 

H-1 

0.40 

Std. 

0 

121.9 

0.5 

H 

H-1 

0.18 

Std. 

0 

121.9 

1.8 

V 

H-1 

0.30 

Std. 

0 

121.9 

1.8 

H 

H-1 

0.14 

Std. 

0 

121.9 

3.7 

V 

H-1 

0.22 

Std. 

0 

121.9 

3.7 

H 

H-1 

0.12 

Std. 

0 

121.9 

6.1 

V 

H-1 

0.15 

Std. 

0 

121.9 

6.1 

H 

H-1 

0.12 

Std. 

120 

33.5 

0.5 

V 

L-1 

3.0 

Std. 

120 

33.5 

0.5 

H 

L-1 

1.8 

Std. 

120 

48.8 

0.5 

V 

L-1 

1.5 

Std. 

120 

48.8 

0.5 

H 

L-1 

0.82 

Std. 

120 

48.8 

1.8 

V 

H-1 

1.1 

Std. 

120 

48.8 

1.8 

H 

H-1 

0.82 

Std. 

120 

48.8 

3.7 

V 

H-1 

0.76 

Std. 
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Table  2.4.  Pre-DICE  THROW  II.  Event  2 (Continued) 
VELOCITY  GAGE  IDENTIFICATION 


Radial 

(deg) 

Range 

(m) 

Depth 

(m) 

Orientation 

Mode* 

Fluid 

Viscosity 

(cs)** 

Gage 

Set 

Range 

(mps) 

Canister 

Type 

120 

48.8 

3.7 

H 

H-1 

0.82 

Std. 

120 

70,1 

0.5 

V 

H-1 

0.76 

Std. 

120 

70.1 

0.5 

H 

H-1 

0.43 

Std. 

120 

70.1 

3.7 

V 

H-1 

0.40 

Std. 

120 

70.1 

3.7 

H 

H-1 

0.43 

Std. 

240 

33.5 

0.5 

V 

L-1 

3.0 

Std. 

240 

33.5 

0.5 

H 

L-1 

1.8 

Std. 

240 

48.8 

0.5 

V 

H-1 

1.5 
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Ranges  in  meters 


Figure  2.2.  Pre-DICE  THROW  II,  Event  2 Gage  Placement  Pattern  Surrounding  Ground  Zero 
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Figure  2.3.  Pre-DICE  THROW  II,  Event  1 Accelerometer  and  Velocity  Gage 
Placement  Profiles 
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(T)  For  Accelerometers  and  Stress  Gages:  DC  op-amp  system;  frequency 
response  flat  to  17  kHz  at  max  gain,  greater  at  lesser  gain. 

For  Velocity  Gages:  3 kHz  carrier-amplifier-demodulator  system; 
frequency  response  flat  to  1 kHz. 

(2)  Digitization  Rates:  6000  samples/sec  for  velocity  and  stress; 
24,000  for  acceleration. 


Figure  2.5.  Pre-DICE  THROW  II  Measurement/Data  Processing  System 
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CHAPTER  3 


PREDICTION  METHODS 

This  chapter  summarizes  the  methods  used  by  the  Air  Force  Weapons 
Laboratory;  Research  and  Development  Associates;  Fie  id  Command,  Defense 
Nuclear  Agency;  and  Waterways  Experiment  Station  in  predicting  ground 
motions  for  the  two  pre-DICE  THROW  II  events. 

A.  AFWL  PREDICTIONS 

The  Air  Force  Weapons  Laboratory  (AFWL)  (Reference  2)  used  the 
predicting  procedures  outlined  in  "Air  Force  Manual  for  Design  and 
Analysis  of  Hardened  Structures,"  AFWL-TR- 74-102,  and  the  results  of 
HE  test  data  analysis  presented  in  "An  Analysis  of  Outrunning  Ground 
Motions,"  AFWL-TR- 74-220. 

The  predicted  crater  volume  was  used  to  predict  the  crater  induced 
motions.  One  departure  from  TR- 74-102  was  made  in  the  displacement 
prediction,  in  that  the  displacement  was  not  assumed  independent  of 
depth.  The  equations  for  near  surface  displacement  were  applied  to 
depths  less  than  0.1  where  V is  predicted  crater  volume,  but  for 

c.  . chs  greater  than  0.1  vV3,  the  displacement  was  attenuated  using 
fits  to  MIDDLE  GUST  data  scaled  to  crater  volume.  Predictions  for 
these  direct  induced  motions  are  given  in  Figures  3.1  through  3.3. 

The  airblast  induced  ground  shock  was  predicted  with  a simplified 
one-dimensional  procedure  recommended  by  TR-74-102  using  a one-dimen- 
sional computer  code.  A bilinear  model  was  assumed  with  Cp  = 0.5 
for  the  soil  above  the  water  table,  and  Cp  = 0.75  C-j  for  the  soil  below 
the  water  table  (Cp  = peak  velocity  and  Ci  = seismic  velocity). 

The  oscillator/  components  of  velocity  at  the  70. 1-,  91. 4-,  and 
121. 9-m  (230-,  300-,  and  400-foot)  ranges  were  predicted  using  the 


^Letter  from  J.  L.  Bratton,  DEVG,  AFWL,  to  Major  T.  Stong,  DNA/SPSS, 
"Empirical  Predictions  for  P*'  *ICE  T.iROW  II,"  12  August  1975. 
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results  of  HE  data  analysis  presented  in  AFWL-TR-74-220.  The  oscil- 
latory components  were  superimposed  with  direct- induced  and  local 
airblast-induced  components  at  these  ranges  according  to  predicted 
time  phasing. 

The  waveforms  were  constructed  from  assumed  airblast  and  direct- 
induced  waveform  components,  which  were  linearly  superimposed.  The 
superposition  time  correlation  was  determined  from  airshock  and  seismic 
propagation  velocities. 

B.  RDA  PREDICTIONS 

R and  D Associates  (RDA)  (Reference  3)  used  an  empirical  approach 
based  on  previous  test  data  to  predict  peak  components  (upward,  down- 
ward, and  outward)  of  velocity  and  displacement  as  functions  of  range 
and  gage  depth.  Comparisons  were  made  of  peak  motion  values  from 
available  surface  tangent  HE  events,  using  yield  to  the  one-third 
power  to  scale  all  distances  and  displacements,  with  velocity  remain- 
ing unsealed. 

The  previous  HE  tests  used  in  predicting  Pre-OICE  THROW  II  ground 
motions  were  MIDDLE  GUST  calibration  shots  1,  4,  5 and  7;  PACE  shots 
1C  and  ID;  MIDDLE  GUST  II,  III  and  IV;  PRAIRIE  FLAT;  DIAL  PACK;  MIXED 
COMPANY  III;  MINE  ORE;  and  DISTANT  PLAIN  6.  Data  were  plotted  from 
these  events  and  curves  fitted  to  these  data  were  then  used  to  construct 
contour  plots  of  peak  velocity  and  displacement  as  functions  of  range 
and  depth.  Figures  3.4  through  3.9  are  the  contour  plots  developed  in 
this  manner.  In  the  case  of  the  velocity  components,  the  fit  was 
generally  mado  through  the  middle  of  the  data.  For  displacements,  it 
was  assumed  that  the  shallow  water  table  would  significantly  enhanc-' 
surface  displacements. 

Estimates  were  made  of  the  periods  associated  with  the  peak  upward, 
downwavd  and  outward  components  of  the  velocity  waveform  using  P = 

Letter  from  Robert  L.  Post,  Jr.,  RDA,  to  Captain  T.  Edwards,  FCONA, 
“Ground  Motion  Predictions  for  Pre-DICE  THROW  II,"  11  August  1975. 
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where  D and  V are  the  peak  displacement  and  velocity  associated  with 
a given  component,  P is  the  period  of  the  velocity  component  and  K 
depends  on  the  shape  of  the  waveform  (K=2  was  used  as  an  upper  bound, 

K=1  for  a lower  bound).  Figure  3.10  gives  the  waveform  predictions 
derived  in  this  manner. 

C.  FIELD  COMMAND  PREDICTIONS 

The  FCDNA  (Reference  4)  predictions  were  made  from  modifications 
of  test  data  recorded  on  the  MIDDLE  GUST  III  event,  which  was  detonated 
in  a geology  not  unlike  the  near  surface  geology  for  the  site  of  the 
Pre-DICE  THROW  II  events.  The  predictions  were  made  with  judgmental 
modifications  to  test  data  resulting  from  differences  in  seismic 
refraction  data  at  the  two  sites.  Consideration  was  also  given  to  the 
differences  in  material  compressibility  at  the  two  sites. 

D.  WES  PREDICTIONS 

Waterways  Experiment  Station  (WES)  (Reference  5)  predictions  were* 
drawn  from  empirical  evaluation  of  data  gathered  from  previous  experi- 
ments at  various  test  sites.  These  experiments  were  MIDDLE  GUST  II, 

III  and  IV  (MG  II,  III  and  IV);  DISTANT  PLAIN  VI  (DP  VI)  and  MINERAL 
ROCK  (MR). 

WES's  predictions  were  made  for  gage  ranging  purposes  only,  and  a 
safety  factor  was  added  so  that  gage  ranges  would  be  adequate  and  the 
signal  conditioning  equipment  could  be  set  with  ample  band  edge  capability. 

The  vertical  and  horizontal  accelerations  were  predicted  using  the 
upper  bound  of  the  peak  data  obtained  from  MG  III,  DP  VI  and  MR  (refer 
to  Figures  3.11  through  3.14).  The  vertical  acceleration  at  the  0.5-m 
(1.5-foot)  depth  followed  the  peak  data  curve  of  MG  III  out  to  the 
36.6-m  (120-foot)  range,  then  followed  the  slope  of  DP  VI  data  out  to 


\etter  from  Captain  T.  Edwards,  FCDNA,  to  Major  T.  Stong,  DNA/SPSS, 
"Pre-DICE  THROW  Predictions,"  28  May  1975. 

^Letter  from  J.  D.  Day,  WES,  to  Major  T.  Stong,  DNA/SPSS,  "Estimates 
of  Ground  Shock  for  Pre-DICE  THROW  II  Series,"  20  May  1975. 
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30^.8  tn  (1000  feet).  At  depths  greater  than  the  0.5  m (1.5  feet), 

MG  III  data  peaks  were  chosen  out  to  39.6  m (130  feet),  then  the 
slope  of  predicted  data  followed  that  of  MR.  The  horizontal  predic- 
tions of  acceleration  at  the  0.5-ni  (1.5-foot)  depth  followed  the 
MG  III  data  out  to  approximately  57.9  m (190  feet),  then  followed  the 
peak  data  obtained  on  the  DP  VI  event.  For  depths  greater  than  0.5  m 
(1.5  feet),  the  MG  III  data  peaks  were  chosen  for  the  prediction. 

Figures  3.15  through  3.17  depict  the  vertical  velocity  predictions 
for  the  0.5-,  1.8-  and  3.7-m  (1.5-,  6-  and  12-foot)  depths,  respec- 
tively. The  shaded  areas  represent  the  data  spread  for  the  crater 
induced  (direct  induced)  motions  from  the  MG  II,  III  and  IV  events. 

Out  to  a range  of  91.4  m (300  feet),  the  point  of  expected  outrunning 
motions,  the  airblast  overpressures  were  expected  to  influence  the 
ground  motions.  A ratio  of  about  0.04  to  0.06  m/sec  per  kpascal  (0.02 
to  0.03  ft/sec  per  psi),  increasing  with  range,  was  used  for  the  pre- 
dictions at  the  0.5-m  (1.5-foot)  depth.  Below  these  depths,  the  pre- 
dictions were  attenuated  to  approximately  70- , 50- , and  30-percent  of 
the  near  surface  velocity  for  the  1.8-,  3.7-,  and  6.1-m  (6-,  12-,  and 
20-foot)  depths,  respectively.  Beyond  the  91.4-m  (300-foot)  range  at 
all  depths,  the  slope  of  the  MG  III  and  IV  upper  boundary  was  followed. 
The  overall  prediction  curve  follows  the  upper  boundary  of  the  MG  data 
with  a 20-  to  30-percent  safety  factor. 

The  horizontal  velocity  predictions  are  shown  in  Figure  3.18. 

There  was  no  attenuation  with  depth  out  to  the  91.4-m  (300-foot)  range. 
Beyond  this  range,  the  velocities  were  attenuated  with  depth  as  shown. 
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Figure  3.2. 
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Figure  3.3.  AFWL  Prediction  of  Vertical  Displacement  versus  Depth  for 
Pre-DICE  THROW  II. 
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Figure  3.4.  RDA  Prediction  of  Airblast  Induced  Vertical  Velocity 
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Figure  3.6.  RDA  Prediction  of  Horizontal  Velocity  for  Pre-DICE  THROW  II 
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Figure  3.8.  RDA  Prediction  of  Direct  Induced  Vertical  Displacement  for  Pre-DICE  THROW  II 


Figure  3.9.  RDA  Prediction  of  Horizontal  Displacement  for  Pre-DICE  THROW  II 
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Figure  3.12.  WES  Prediction  of  Vertical  Acceleration  at  Depths 
Greater  than  0.5  m for  Pre-DICE  THROW  II. 
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CHAPTER  4 


COMPARISON  OF  DATA  AND  PREDICTIONS;  EVENT  1 


The  purpose  of  this  chapter  is  to  compare  the  predictions  des- 
cribed in  Chapter  3 with  the  actual  data  obtained  from  Event  1.  Simi- 
lar comparisons  for  Event  2 will  be  made  in  Chapter  5.  An  effort  has 
been  made  to  present  the  comparisons  in  a way  that  will  be  as  meaning- 
ful as  possible.  It  is  hoped  that  they  will  prove  useful  as  a reference 
for  future  HE  test  predictions  and  data  analysis. 


A.  TABULATIONS  OF  DATA 


All  of  the  numerical  peak  values  used  in  this  series  of  comparisons, 
both  experimental  and  predicted,  are  given  in  Tables  4.1  through  4.3. 

Tables  4.1.  4-2,  and  4.3  contain  the  acceleration,  velocity,  and  dis- 
placement V •,  respectively,  for  Event  1.  If  more  than  one  measure- 
ment was  maae  at  the  same  radial  distance  from  ground  zero  and  the  same 
depth,  but  on  a different  radial,  then  the  average  value  of  these  data 
is  computed,  and  these  averages  are  also  given  in  the  tables.  Additional 
details  pertaining  to  transducer  type  and  placement  are  given  in  Chapter  2. 
The  displacement  data  given  in  Table  4.3  are  integrations  of  the  veloc- 
ity records,  and  are  included  only  to  indicate  trends  relative  to  the 
predictions. 

It  should  be  noted  that,  although  the  data  columns  in  the  tables 
are  titled  ''Direct  Induced  (Positive)"  and  “Airblast  Induced  (negative)", 
these  designations  are  correct  only  for  the  vertically  oriented  gages. 

In  the  case  of  the  horizontal  gages,  peak  positive  excursions  of  the 
traces  were  recorded  in  the  "Direct  Induced"  column,  and  peak  negative 
excursions  were  recorded  in  the  "Airblast  Induced"  column.  No  attempt 
was  made  to  separate  the  two  types  of  motions  in  the  horizontal  data, 
and  the  recorded  negative  peaks  probably  arise  from  elastic  recovery 
of  the  medium  in  which  the  horizontal  gage  is  situated. 

In  the  prediction  columns,  the  peak  values  predicted  by  each  agency 
for  direct  induced  (DI)  or  positive  excursions  and  for  airblast  induced 
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(AI)  or  negative  excursions  are  given.  Where  a prediction  was  not 
available  at  the  specified  location,  the  symbol  or  is  used. 

Deviations  of  the  predictions  from  the  experimental  data  averages 
were  also  computed.  The  percentage  deviations  were  determined  from  the 
following  formula: 


Percent  Deviation  (n) 


Average  Experimental  Value  - Prediction 
Average  Experimental  Value  ^ 


A positive  deviation  indicates  that  the  average  data  point  was  greater 
than  that  predicted,  whereas  a negative  deviation  indicates  that  it 
was  less.  It  is  obvious  that  the  maximum  positive  percentage  deviation 
is  100,  while  there  is  no  limit  to  the  negative  deviation. 

Only  three  types  of  data  are  included  in  this  report,  acceleration 
and  velocity  measurements  and  the  integrated  velocity  data,  or  displace- 
ment. Some  predictions  of  stress  were  made  by  FCDNA  and  AFWL,  but 
these  predictions  were  for  ranges  and  depths  not  covered  by  the  trans- 
ducer array.  Hence  no  attempt  has  been  made  to  ■incorporate  considera- 
tions of  stress  into  the  present  study. 

B.  GRAPHICAL  PRESENTATIONS 

Following  the  tables,  a large  number  of  plots  are  given  showing 
the  experimental  data  and  predictions  in  different  forms.  The  first 
group  of  plots.  Figures  4.1  through  4.6,  shows  the  vertical  accelera- 
tion data  versus  range  from  ground  zero.  There  is  a separate  plot  for 
each  depth  at  which  experimental  data  and  predictions  are  available. 

The  direct  induced  peaks  are  plotted  as  positive  values,  whereas  the 
airblast  induced  peaks  are  plotted  as  negative  values.  Calculated 
percentage  deviations  of  the  predictions  from  the  data  are  also  plotted 
separately  for  the  direct  and  airblast  induced  cases.  In  Figures  4.7 
through  4.9,  the  same  information  is  plotted  as  functions  of  depth,  with 
a separate  plot  for  each  range. 

Horizontal  accelerations  versus  range  are  given  in  Figures  4.10 
through  4.12.  The  same  data  versus  depth  is  plotted  in  Figures  4.13 
through  4.15.  It  should  be  noted  that,  in  these  plots  as  well  as  in 
others  in  this  and  the  following  chapter,  only  those  data  points  for 
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which  a comparison  between  prediction  and  experiment  was  possible  have 
beer  included.  Predictions  at  locations  for  which  no  experimental  data 
was  obtained  have  been  omitted,  as  have  experimental  points  at  locations 
for  which  no  predictions  were  made.  In  a few  cases,  complete  graphs 
have  been  omitted  when  there  were  no  comparisons  to  be  made. 

Beginning  with  Figure  4.16,  the  velocity  comparisons  are  given. 

Where  an  average  velocity  value  was  computed  because  there  was  more  than 
one  gage  at  the  designated  range  and  depth,  this  average  value  is  plotted 
with  a vertical  bar  indicating  the  actual  data  spread.  A single  point 
with  no  bar  means  either  that  only  one  value  was  obtained  or  that  the 
data  spread  was  not  detectable  on  the  scale  of  the  plots.  Actual  values 
may  be  obtained  by  referring  to  the  tabulations.  Figures  4.16  through 
4.21  present  vertical  velocities  versus  range  for  each  depth,  together 
with  deviations  of  the  predictions,  and  Figures  4.22  through  4.25  pre- 
sent vertical  velocities  versus  depth  for  each  range.  Horizontal  veloc- 
ity peak  values  versus  range,  with  deviations  of  the  predictions,  are 
given  in  Figures  4.26  through  4.28,  and  the  same  values  versus  depth 
are  given  in  Figures  4.29  through  4.32. 

Finally,  the  displacement  peak  values  are  given  in  Figures  4.33 
through  4.48.  As  mentioned  previously,  experimental  values  are  taken 
from  integrations  of  the  velocity  traces,  so  the  data  shown  here  is  not 
independently  obtained.  Again,  averages  of  the  data  points  are  used 
where  there  was  more  than  one,  with  vertical  bars  to  indicate  the  data 
spread.  Vertical  displacements  versus  range,  with  separate  prediction 
deviation  plots,  are  shown  in  Figures  4.33  through  4.37;  vertical  dis- 
placements versus  depth  are  shown  in  Figures  4.38  through  4.41.  The 
horizontal  displacement  data  is  given  as  functions  of  range  in  Figures 
4.42  through  4.44,  and  as  functions  of  depth  in  Figures  4.45  through  4.48. 

C.  WAVEFORM  COMPARISONS 

In  a few  cases,  attempts  were  made  by  Field  Command  (FC)  and  the 
Air  Force  Weapons  Laboratory  (AFWL)  to  predict  the  complete  waveform 
traces  to  be  expected  from  the  gages.  This  involves  superposition  of 
the  expected  direct  induced  and  airblast  induced  motions,  with  appropriate 


pulse  durations  of  each.  The  phase  relationship  of  the  two  signals, 
or  time  of  arrival  of  each,  must  also  be  predicted. 

These  predicted  waveforms  are  compared  with  actual  recorded  traces 
from  Event  1,  where  both  are  available,  in  Figures  4.49  through  4.62. 
Although  some  features  of  the  waveforms  are  predicted  with  remarkable 
accuracy,  it  is  obvious  that  our  understanding  of  all  aspects  of  ground 
motion  is  not  sufficient  that  reliable  waveform  predictions  can  be  made 
in  every  case. 

D.  OBSERVATIONS  ON  EVENT  1 COMPARISONS 

In  Tables  4.4,  4.5,  4.6,  and  4.7,  observations  from  the  comparisons 
of  data  and  predictions  shown  in  the  plots  are  made  separately  for  the 
four  predicting  agencies.  In  some  cases  the  amount  of  comparable  data 
was  too  sparse  to  indicate  any  obvious  trends.  These  cases  are  noted 
in  the  tables.  A final  observation  summary  covering  all  agencies  is 
given  in  Table  4.8. 

It  should  be  mentioned  here  that  there  is  no  intent  in  this  summary 
to  imply  that  any  one  set  of  predictions  is  better  or  worse  than  any 
other.  It  is  recognized  that  the  predictions  were  made  for  different 
purposes  (e.g.,  for  gage  range  setting)  and  with  varying  degrees  of 
effort  involved.  Rather,  it  is  hoped  that  the  reader  will  gain  an  over- 
all feeling  for  the  accuracy  with  which  ground  motion  predictions  can 
be  made,  and  an  appreciation  for  the  care  that  must  be  taken  to  prepare 
a good  prediction.  This  feeling  and  appreciation  should  serve  as  back- 
ground for  the  preparation  and  evaluation  of  predictions  for  future 
tests. 

1.  AFWL  Predictions.  The  Air  Force  Weapons  Laboratory  predicted 
airblast  induced  vertical  and  horizontal  accelerations,  velocities,  and 
displacements.  They  also  predicted  direct  induced  vertical  velocities 
and  displacements.  Observations  on  comparisons  of  these  predictions 
with  experimental  data,  along  with  references  to  the  figures  from 
which  the  observations  were  made,  are  given  in  Table  4.4.  The  AFWL 
predictions  for  direct  induced  motions  were  generally  higher  than  the 
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data.  It  will  be  recalled  that  these  predictions  were  based  on  an 
estimate  of  anticipated  crater  volume.  The  predicted  crater  volume 
was  105,000  cubic  feet,  whereas  actual  crater  volume  was  152,000  cubic 
feet.  It  may  be  speculated  that  the  AFWL  predictions  would  have  been 
even  higher,  and  hence  further  from  the  measured  values,  if  the  crater 
volume  prediction  had  been  closer. 

2.  RDA  Predictions.  R and  D Associates  predicted  airblast  in- 
duced and  direct  induced  vertical  velocities  and  displacements.  Pre- 
dictions were  also  made  of  horizontal  velocity  and  displacement  peaks. 
These  latter  were  not  specified  as  arising  from  airblast  or  direct 
induced  motions.  Observations  on  comparisons  of  these  predictions 

with  the  experimental  data  are  given  in  Table  4.5.  The  airblast  induced 
motions,  as  predicted  by  RDA,  did  not  take  into  account  the  influence 
of  the  water  table  at  a depth  of  approximately  two  meters  (6-8  feet). 
The  experimental  data,  however,  show  a substantial  change  in  airblast 
induced  motions  across  this  boundary,  and  it  shows  up  on  the  data  and 
deviation  plots.  The  RDA  predictions  are  generally  low  above  the  water 
table  and  high  below  it. 

3.  FC  Predictions.  Field  Command  predicted  only  direct  induced 
vertical  and  horizontal  velocities.  Very  few  predictions  were  made, 
and  these  were  intended  only  as  spot  checks  as  to  whether  or  not  other 
predictions  were  reasonable.  The  approach  was  empirical  and  intuitive; 
however,  results  were  excellent  in  that  they  came  quite  close  in  most 
cases  to  the  measured  velocity  peaks.  Observations  on  the  comparisons 
are  given  in  Table  4.6. 

4.  WES  Predictions.  Waterways  Experiment  Station  predicted 
direct  induced  vertical  and  horizontal  accelerations  and  velocities. 
Observations  on  comparisons  of  these  predictions  with  the  experimental 
data  appear  in  Table  4.7.  It  would  be  expected  that  these  predictions, 
because  they  were  made  for  the  purpose  of  setting  the  gage  ranges, 
would  tend  to  be  conservative;  that  is,  somewhat  higher  than  the  meas- 
ured peak  values.  Study  of  the  plotted  comparisons  indicates  chat  this 
is  indeed  the  case. 


Table  4.1.  Pre-DICE  THROW  II,  Event  1 Acceleration  Data  with  Prediction  and  Deviations  (Continued) 


Table  4.1.  Pre-DICE  THROW  II,  Event  1 Acceleration  Data  with  Prediction  and  Deviations  (Continued) 


Table  4.1.  Pre-DICE  THROW  II,  Event  1 Acceleration  Data  with  Preoiction  and  Deviations  (Continued) 


Table  4.2.  Pre-DICE  THROW  II,  Event  1 Velocity  Data  with  Predictions  and  Deviations  (Continued) 


^ eg  lA 

r-  VO  0«  lO 

III  I 


1 

in 

m 

o 

1 

en 

OD 

f 

IN. 

lO 

o 

t 

o 

o 

1 

•Vk 

'V.. 

■Nw 

•N.. 

o\ 

a-. 

V£> 

CO 

in 

m 

m 

ifi 

m 

o 

o 

o 

o 

O 

o 

o 

^i**^*^  ooo  egoeg  OO 


ooo  oovo  o o e 
o 


,~y  W I ^ «l 

4»?i+>>(l)  COO>CT» 

t.  .2  -r-  ... 


.2'^ 

■o  iy  o 

n3  3 
a:  ^ 


S 5 o 

r—  OJ 


CNIr-vO  lOrOO  eg.—  r«-  mO^CNJ  r-  . If) 

ooo^  egcg«^  ^egr^.  ^vo*cf  rgroro  ro^ro  O ^ 

• • ‘ ...  . 

o ^ OOO  ooo  OOO  OOO  ooo  o •—  o 


OS  OS  C>s  OO  OQ  OO 

ocsi^  ocgv  oeg^  opg^  oeg*^ 

.—  eg  .—  eg  .—  eg  — eg  ^eg  .—  eg 


a:xx  xxx 


Table  4.2.  Pre-DICE  THROW  II,  Event  1 Velocity  Data  with  Predictions  and  Deviations  (Continued) 
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Table  4.2.  Pre-DICE  THROW  II,  Event  1 Velocity  Data  with  Predictions  and  Deviations  (Continued) 
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Table  4.3.  Pre-DICE  THROW  II,  Event  1 Displacement  Data  with  Predictions  and  Deviations 


Acceleration  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  1 
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re  4.12.  Horizontal  Acceleration  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  1 
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4.16.  Vertical  Velocity  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  1 
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Vertical  Velocity  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  1 


igure  4.18.  Vertical  Velocity  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  1 
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jre  4.20.  Vertical  Velocity  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  1 
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4.21.  Vertical  Velocity  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  1 
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Figure  4,27.  Horizontal  Velocity  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  I 


Figure  4.28.  Horizontal  Velocity  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  1 


Figure  4.33,  Vertical  Displacement  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  1 
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Figure  4.38.  Vertical  Displacements  vs.  Depth  with  Predictions  for  Pre-DICE  THROW  II,  Event  1 
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4.42.  Horizontal  Displacement  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  1 


Figure  4.43.  Horizontal  Displacement  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  1 
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Horizontal  Displacement  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  1 


Figure  4.46.  Horizontal  Displacement  vs.  Depth  with  Predictions  for  Pre-DICE  THROW  II,  Event  1 
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Figure  4.47.  Horizontal  Displacement  vs.  Depth  with  Predictions  for  Pre-DICE  THROW  II,  Event  1 


-2* 


0°  RADIAL 


ACTUAL 

FC 


EVENT  1 

VERTICAL  VELOCITY 
24.2  m RANGE 
.5  m DEPTH 


120°  RADIAL 


240°  RADIAL 


X7 


TIME  FROM  DETONATION  (sec) 


Figure  4.49.  Velocity  Waveforms  and  Predictions  for  Pre-DICE  THROW  II,  Event  1 
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Figure  4.60.  Velocity  Waveforms  and  Predictions  for  Pre-DICE  THROW  II,  Event  1 


VERTICAL  VELOCITY  (m/sec) 


ACTU/ 

0°  RADIAL  AFWL 


.5  m DEPTH 


120°  RACIAL 


HORIZONTAL  VELOCITY  (m/sec) 


4-83 


As  compared  to  data  from  Event  1,  the  FC/DNA  predictions  were. 


CHAPTER  5 


COMPARISON  OF  DATA  AND  PREDICTIONS;  EVENT  2 

This  chapter  is  identical  to  Chapter  4 except  that  data  from  the 
120-ton  ANFO  detonation,  Event  2,  was  used  as  a base  against  which 
comparisons  were  made.  The  predictions  in  all  cases  were  the  same 
for  both  events,  so  any  differences  between  the  comparisons  presented 
here  and  those  of  Chapter  4 will  be  entirely  due  to  differences  in  the 
experimental  data  base. 

A.  TABULATED  AND  GRAPHICAL  PRESENTATIONS 

The  material  here  is  organized  exactly  the  same  as  that  in  Chapter 
4,  and  many  of  the  comments  made  in  that  chapter  also  apply  here. 
Therefore,  the  reader  is  advised  to  look  at  Chapter  4 before  proceeding. 

Experimental  and  predicted  peak  values  for  acceleration,  velocity, 
and  displacement,  respectively,  are  given  in  Tables  5.1,  5.2,  and  5.3. 
Percentage  deviations,  calculated  as  previously  described,  are  also 
given  in  the  tables.  As  there  were  fewer  ground  motion  data  recorded 
on  Event  2 than  on  Event  1,  the  tables  are  somewhat  shorter  than  those 
in  Chapter  4. 

Graphical  presentations  of  the  tabulated  data  comprise  Figures  5.1 
through  5.39.  As  before,  in  the  case  of  vertical  motions,  an  effort 
is  made  to  separate  the  direct  induced  (positive)  and  airblast  induced 
(negative)  motions.  The  horizontally  oriented  gages  record  both  types 
of  motion  in  the  positive  (outward)  direction,  and  in  these  cases  no 
attempt  is  made  to  separate  the  two  motions.  Peak  positive  and  negative 
excursions  of  the  traces  are  'recorded. 

Vertical  acceleration  data  and  predictions  are  plotted  versus  range 
in  Figures  5.1  through  5.4.  Deviations  of  the  predictions  are  also 
given.  The  same  data  are  given  as  functions  of  depth  in  Figures  5.5 
and  5,6.  The  corresponding  horizontal  acceleration  data  is  plotted 
against  range  in  Figures  5.7  through  5.9  and  against  depth  in  Figures 
5.10  and  5.11. 
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Figures  5.12  through  5.16  present  the  vertical  velocity  data  and 
deviations  for  Event  2 versus  range.  . igures  5.17  through  5.19  show 
this  same  data  as  functions  of  depth.  The  horizontal  velocity  data  is 
given  versus  range  in  Figures  5.20  through  5.22  and  versus  depth  in 
Figures  5.23  through  5.25. 

Displacement  predictions  and  integrations  of  the  velocity  traces 
are  also  presented  graphically,  as  are  the  calculated  deviations  for 
this  data.  Vertical  displacements  are  plotted  versus  range  in  Figures 
5.26  through  5.30,  and  versus  depth  in  Figures  5.31  through  5.33. 
Horizontal  displacements  versus  range  are  presented  in  Figures  5.34 
through  5.36.  The  same  horizontal  displacements  appear  as  functions 
of  depth  in  Figures  5.37  through  5.39. 

A few  comparisons  of  predicted  and  experimentally  obtained  wave- 
forms are  also  possible  for  the  data  of  Event  2.  These  comparisons  are 
restricted  to  vertical  and  horizontal  velocity  traces.  As  before,  cer- 
tain features  of  the  records  are  remarkably  well  predicted,  while  others 
leave  much  to  be  desired.  The  traces  are  reproduced  in  Figures  5.40 
through  5.52. 

B.  OBSERVATIONS  ON  EVENT  2 COMPARISONS 

Tables  5.4  through  5.7  are  used  to  present  observations  on  the 
graphical  representations  of  the  data  comparing  predictions  with  Event  2 
data.  The  format  is  exactly  as  that  of  Chapter  4,  and  comments  made  in 
that  chapter  also  apply  to  the  present  comparisons.  Finally,  a summary 
of  the  observations  covering  predictions  made  by  all  four  agencies  is 
given  in  Table  5.8.  It  will  be  seen  from  comparison  of  Table  5.8  with 
Table  4.8  that  the  relationship  of  predictions  to  the  data  was  almost 
the  same  for  both  events.  The  actual  crater  volume  for  this  event  was 
166,000  cubic  feet,  as  compared  to  152,000  cubic  feet  for  Event  1 and 
to  105,000  cubic  feet  predicted  by  AFWL. 
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Table  5.2.  Pre-DICE  THROW  II,  Event  2 Velocity  Data  with  Predictions  and  Deviations  (Continued) 
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Table  5.3.  Pre-DICE  THROW  II,  Event  2 Displacement  Data  with  Predictions  and  Deviations  (Continued) 


Table  5.3.  Pre-DICE  THROW  II,  Event  2 Displacement  Data  with  Predictions  and  Deviations  (Continued) 
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Figure  5.3.  Vertical  Acceleration  vs.  Range  with;  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event 
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Figure  5.5.  Vertical  Acceleration  vs.  Depth  with  Predictions  for  Pre-DICE  THROW  II,  Event  2 


PREDICTION  LESS 
THAN  ACTUAL 


PftEDICTIOM  GREATER 
THAN  TCTUAL 


Range,  meters 

Figure  5.12.  Vertical  Velocity  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  2 
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Figure  5.13.  Vertical  Velocity  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  2 
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Figure  5.14,  Vertical  Velocity  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  2 
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Figure  5.15.  Vertical  Velocity  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  2 


Vertical  Velocity  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  2 


Figure  5.17.  Vertical  Velocity  vs.  Depth  with  Predictions  for  Pre-DICE  THROW  II,  Even 


Figure  5.19.  Vertical  Velocity  vs.  Depth  with  Predictions  for  Pre-DICE  THROW  II,  Event  2 


Figure  5.23.  Horizontal  Velocity  V5.  Depth  with  Predictions  for  Pre-DICE  THROW  II,  Event  2 
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Vertical  Displacement  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  2 


Figure  5.27.  Vertical  Displacement  vs.  Range  with  Predictions  and  Deviations  for  Pre-DICE  THROW  II,  Event  2 
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5.29.  Vertical  Displacement  vs.  Range  with  Predictions  and  Deviations  foi-  Pre-DICE  THROW  II,  Event  2 


Vertical  Displacement  vs.  Depth  with  Predictions  for  Pre-DICE  THROW  II,  Event  2 


Figure  5.33.  Vertical  Displacement  vs.  Depth  with  Predictions  for  Pre-DICE  THROW  II,  Event  2 
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Figure  5.49.  Velocity  Waveforms  and  Predictions  for  Pre-DICE  THROW  II,  Event  2 
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Table  5.5.  Observations  on  Comparisons  of  P'  ’’redictions  with  Event  2 Data 


As  compared  to  data  from  Event  2,  the  FC/ONA  predictions  were: 


Table  5.7.  Observations  on  Conparisons  of  WES  Predictions  with  Event  2 Data 


Table  5.8.  Data  versus  Prediction  Summary  for  Pre-DICE  THROW  II,  Event  2 


CHAPTER  6 


COMPARISON  OF  EVENTS  1 AND  2 DATA 

A small  effort  has  been  made  to  compare  the  experimental  results 
of  Event  1 with  those  of  Event  2.  Although  much  larger  effort  is  jus- 
tified in  order  to  understand  cratering  phenomena,  the  scope  of  this 
analysis  does  not  permit  extensive  consideration.  Consequently,  the 
present  effort  looks  only  at  those  aspects  of  comparison  between  the 
two  events  which  might  bear  on  predictive  capabilities. 

It  might  be  expected  that,  because  the  charge  yields  were  similar 
(however,  the  energy  releases  are  not  yet  known),  the  far-field  effects 
would  also  be  similar.  In  that  the  charge  configurations  were  dif- 
ferent and  the  fireball  photos  revealed  a somewhat  different  early- 
time  development  in  the  two  cases,  different  near- field  effects  may 
be  expected.  The  craters  were  indeed  measurably  different. 

The  water  table  roughly  defined  the  interface  between  low  and 
high  wave  propagation  velocity  layers,  as  described  in  Chapter  1. 

Because  the  crater  of  Event  1 penetrated  the  water  table,  it  is  clear 
that  a significant  part  of  the  blast  energy  was  injected  into  the  lower 
layer.  The  reverse  was  true  for  Event  2,  in  which  less  energy  was  in- 
jected into  the  lower  layer  and  the  crater  was  broader  and  shallower. 
These  comments  are  qualitatively  supported  by  observations  from  the 
data  as  indicated  in  Table  6.1. 

More  detail  is  available  in  the  illustrations  of  Figures  6.1 
through  6.17.  Although  the  data  is  scant,  it  can  be  seen  that  measure- 
ments are  in  most  cases  very  similar  for  the  two  shots.  In  a few 
cases,  notably  the  air- induced  vertical  velocities  at  1.8  m depth  (che 
negative  peaks  shown  in  the  lower  half  of  Figure  6.6),  the  discrep- 
ancies in  measurement  are  large.  A comparable  effect  at  the  same 
ranges  at  the  shallower  depth  (O.B  meter)  is  not  observed.  Because 
Events  1 and  2 were  not  at  the  same  site,  this  anomaly  could  be  one 
of  geology.  The  horizontal  velocity  data  at  these  same  stations 
(Figure  6.9),  however,  shows  far  less  divergence  between  Events  1 and  2. 
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It  Is  possible  that  an  instrumental  problem  could  be  responsible  for 
the  differences. 

In  summary,  the  differences  detectable  in  ground  motion  measure- 
ments far  from  the  craters  of  Events  1 and  2 are  minor,  except  in  a 
few  cases.  Measurements  in  the  vicinity  of  the  craters  should  show 
diffe»^ences  due  to  charge  shape.  Comparison  of  Figures  4.17  with 
5.13  and  4.34  with  5.27  illustrates  a slightly  better  predictive 
capability  for  the  tangent  sphere  (Event  1)  charge  configuration, 
although  there  is  not  enough  data  to  draw  a strong  conclusion.  In 
these  cases,  predictions  (by  RDA)  were  based  on  scaling  data  from 
tests  which  were  of  the  tangent  sphere  configurations.  Better  predic 
tive  methods,  or  a more  extensive  data  base  for  configurations  other 
than  the  tangent  sphere,  may  be  needed  to  distinguish  charge  shape 
effects. 
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Table  6.1.  Observations  on  Comparisons  of  Data  from  Events  1 and  2 
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Figure  6.9.  Horizontal  Velocity  vs.  Range  for  Pre-DICE  THROW, 
Events  1 and  2 
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CHAPTER  7 


ANALYSIS  AND  CONCLUSIONS 


A.  INTRODUCTION 

The  objective  of  this  report  was  to  compa'^e  and  to  evaluate,  if 
possible,  the  several  prediction  techniques.  Casual  inspection  of  the 
results  does  not  provide  a clear-cut  answer  to  the  question  of  which 
technique  most  accurately  predicted  the  results  obtained.  A more 
quantitative  approach,  using  some  simple  statistical  methods,  is  at- 
tempted in  this  chapter.  The  data  base  does  not  really  support  a sta- 
tistical analysis,  but  casting  the  numbers  in  statistical  terminology 
does  provide  a familiar  and  useful  approach  to  comparisons. 

In  making  the  analysis,  the  particle  velocity  data  are  used  be- 
cause these  measureti/^nts  are  numerous,  are  of  good  quality,  and  are 
relatively  easy  to  interpret.  There  is,  however,  a limitation  which 
affects  the  validity  of  the  present  discussion.  In  reading  the  data 
traces,  the  main  positive  and  negative  peaks  (one  each)  were  read. 

For  the  vertical  velocity  components,  these  peaks  invariably  occurred 
within  the  first  100  milliseconds  of  the  trace,  and  corresponded 
respectively  to  the  leading  edges  of  the  air-  and  direct- induced  motions. 
For  the  horizontal  velocity  components,  however,  the  situation  is  not 
so  easily  interpreted.  The  inclination  of  the  air-  and  direct- induced 
waves  close-in  to  the  burst  is  very  large  and  the  horizontal  or  radial 
components  of  the  particle  motion  are  therefore  relatively  small.  Also, 
both  air-  and  direct- induced  motions  are  in  the  outward,  or  positive, 
direction.  Consequently  the  peaks  read  were  frequently  due  not  to 
air-  and  direct-induced  motions,  but  to  another  type  of  wave,  the 
surface  wave  (probably  of  the  Raleigh  type)  which  results  from  sudden 
depression  of  the  ground  at  the  burst  point.  These  peaKS  invariably 
occur  farther  out  in  time,  typically  at  200  - 300  milliseconds.  Seen 
on  the  truces,  the  slow  undulations  of  the  surface  waves  are  easy  to 
recognize.  In  the  tables,  it  cannot  be  determined  whether  a given 
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peak  value  arose  from  airblast-induced,  direct-induced,  or  surface 
wave  motions. 

The  long  duration  surface  wave  motions  were  more  significant  in 
these  events  than  had  been  anticipated  for  data  analysis.  Beyond  a 
range  of  about  50  meters,  however,  there  is  a tendency  for  the  peak 
values  of  the  induced  waves  to  dominate  the  peaks  of  the  surface  waves. 
Nevertheless,  it  is  still  not  possible  to  distinguish  entries  in  the 
tables  according  to  whether  they  are  air-  or  direct-induced.  Conse- 
quently, we  believe  it  is  impractical  to  make  prediction  comparisons 
on  the  basis  of  horizontal  peak  values  until  the  data  is  rescanned  under 
the  selective  intent  to  pick  out  the  airblast-  and  direct-induced  peaks. 

The  vertical  velocity  peaks,  as  previously  explained,  can  almost 
always  be  interpreted  unambiguously.  There  are  a few  cases  where  iden- 
tification of  the  wave  could  be  in  doubt,  but  these  cases  will  prob- 
ably not  comprise  a significant  fraction  of  the  total.  Hence  in  the 
present  analysis,  the  vertical  velocity  peaks  will  be  used  as  given. 

Most  of  the  measurements  on  Events  1 and  2 were  taken  in  the  low 
velocity  layer  where  the  "X"  discussed  in  Chapter  1 is  found.  Conse- 
quently they  are  typical  of  regions  where  both  air-  and  direct- induced 
waves  are  important.  If  the  first  peak  is  negative,  the  air-induced 
wave  arrived  first,  and  if  the  first  peak  is  positive,  then  the  direct- 
induced  wave  arrived  first.  One  of  the  important  tasks  of  the  predic- 
tion efforts  was  to  try  to  estimate  the  contributions  of  these  two 
waves  individually.  Therefore,  the  vertical  velocity  data  are  ideal 
suited  to  the  analysis  undertaken  here.  It  will  be  obvious  that  the 
analysis  is  far  from  exhaustive,  and  that  a re-reading  of  the  data 
traces,  together  with  a more  sophisticated  analytical  approach,  would 
undoubtedly  achieve  more  authoritative  results. 

B . SIMPLE  AVERAGES  OF  PERCENT  DEVIATIONS 

We  look  first  at  simp'^e  averages  of  the  percent  deviations  appear- 
ing in  Chapters  4 and  5 to  try  to  find  trends  in  the  easiest  way.  The 
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data  to  be  used  are  the  observed  peak  values  compared  to  the  predicted 
peak  values,  called  "percent  deviations"  and  computed  from 


Pe'^cent  Deviation  (n) 


Average  Experimental  Value  - Prediction 
Average  Experimental  Value  ^ 


'b  a statistical  variable,  n is  unusual  because  its  positive  range 
^ truncated;  that  is,  all  positive  values  must  fall  between  0 and  100 
while  the  negative  values  are  distributed  between  0 and  -oo.  Conse- 
quently using  these  data  leads  to  somewhat  skewed  averages,  especially 
when  there  are  numerous  values  in  the  sample  of  magnitude  less  than 
-100.  The  small  sample  does  not  justify  a more  careful  analysis,  how- 
ever. For  our  purposes,  we  believe  simple  comparisons  are  sufficient 
to  draw  the  conclusions  called  for. 

The  simple  averages  of  the  percent  deviations  are  computed  in  the 
following  fashion: 


where  the  are  the  percent  deviations  at  all  depths  for  a given  range, 
and  the  sign  of  is  taken  into  account.  The  results  are  given  in 
Tables  7.1  and  7.2.  The  Field  Command  predictions  were  not  used  in  this 
analysis  because  the  number  of  available  predictions  was  so  small. 

One  sees,  first,  a preponderance  of  negative  values,  which  signi- 
fies a tendency  to  overestimate  the  magnitude  of  the  motions.  This 
tendency  seems  to  increase  with  range.  An  exception  is  RDA's  predic- 
tions of  the  DI  peak  values,  for  which  the  averages  of  the  deviations 
are  nearly  all  positive.  The  stated  purpose  of  the  WES  predictions, 
as  mentioned  earlier,  was  for  gage  ranging,  whicn  requires  a safe  over- 
estimate in  order  to  avoid  loss  of  data.  Therefore,  the  WES  averages 
are  all  negative,  as  they  should  be. 

Second,  one  observes  that  magnitudes  of  almost  all  the  n are 
greater  than  50%.  In  other  words,  the  probability  of  predicting  ground 


Table  7.1.  Averages  of  Percent  Deviations  at  Constant 
Ranges  for  Vertical  Velocity  Predictions; 
Event  1 


Range  (m) 

AFWL 

WES 

RDA 

DI 

AI 

DI 

DI 

AI 

24.4 

-36.1 

-12.1 

-46 

27.5 

28.7 

33.5 

- 

- 

-48.5 

29.6 

-39.5 

48.8 

-280 

-134 

-65.8 

3.1 

-143 

70.1 

-57.4 

-268 

-58.3 

3.0 

-161 

91.4 

-91.5 

-126 

-75.0 

1.4 

-126 

121.9 

-89  5 

-73.9 

-109 

6.4 

-198 

182.9 

- 

- 

-415 

-67.4 

-44.5 

Table  7.2.  Averages  of  Percent  Deviations  at  Constant 
Ranges  for  Vertical  Velocity  Predictions; 
Event  2 


Range  (m) 

AFWL 

WES 

RDA 

DI 

AI 

DI 

DI 

AI 

24.4 

29.1 

-9.2 

-34.3 

32.3 

-4.9 

33.5 

- 

- 

19.1 

-4.4 

48.8 

-35.8 

-139 

8.5 

-28.5 

70.1 

-91.5 

-80.2 

-11.2 

-50.3 

91.4 

- 

-5.5 

-22.5 

42.5 

44.5 

121.9 

-70.5 

-15.4 

-135 

-10.8 

-94.1 

motion  peaks  accurately  is,  on  the  average,  no  better  than  50%.  A 
more  sensitive  measure  of  this  quantity  is,  however,  shown  in  the  next 
section. 

C.  STANDARD  DEVIATIONS 

The  above  simple  avp>'V!ges  do  not  really  provide  enough  insight 
into  the  acci'racy  of  the  predictions  so  we  now  examine  standard  devia- 
tions. The  statistical  deviation  of  the  members  of  a set  from  the 
centroid  of  the  set  is  an  approved  measure  of  the  scatter,  or  of  the 
accuracy  in  this  case,  of  the  members  of  the  set.  We  again  calculate 
simple  averages  of  the  percent  deviations  appearing  in  Chapters  4 and 
5 from  the  expression: 

n 


This  time,  however,  the  set  of  are  chosen  for  all  ranges  for  a given 
depth.  This  choice  results  from  the  conviction  that  the  presence  of 
the  low  wave  velocity  surface  layer  plays  a dominant  role  in  determin- 
ing the  observed  ground  motion  phenomenology.  The  n computed  here  will 
not  be  the  same  as  the  n used  in  Section  B because  the  sets  over  which 
the  averages  are  taken  are  different. 

Standard  deviations  of  the  set  of  percent  deviations  from  these 
average  values  can  then  be  calculated  from  the  usual  formula: 


a 


1/2 


As  before,  the  values  of  n are  not  uniformly  distributed  about  zero; 
positive  values  cannot  be  larger  than  100.  The  results  in  a will  also 
be  skewed.  This  will  not  be  significant,  however,  for  the  intuitive 
interpretations  to  be  drawn  here. 

The  averages  of  the  percent  deviations  at  constant  depths,  and 
their  standard  deviations,  are  given  in  Table  7.3.  This  presentation 
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of  the  data  allows  an  evaluation  of  how  the  three  prediction  techniques 
account  for  the  layer  effect.  Each  entry  in  the  table  consists  of 
two  terms,  in  the  form  n ± a. 


Table  7.3.  Averages  of  Percent  Deviations  at  Constant 
Depths  with  Standard  Deviations 


Event  1 

AFWL 

WES 

RDA 

Depth  (m) 

DI 

AI 

DI 

DI 

AI 

0.5 

-0.94±28 

-45±44 

-130±155 

0.53153 

16131 

1.8 

- 

- 

-1751106 

-32134 

-191107 

3.7 

-179±63 

-170±121 

-1291110 

-5127 

-2101127 

6.1 

- 

-205±141 

-31121 

2218 

-196187 

9.1 

- 

- 

3713 

62113 

-3421285 

Event 

? 

AFWL 

WES 

RDA 

Depth  (m) 

DI 

AI 

DI 

DI 

AI 

0.5 

6.9±33 

15±14 

-64180 

33116 

1.8 

- 

- 

-124153 

-35191 

3.7 

-177±62 

-60±49 

-102117 

-10119 

-107143 

6.1 

- 

-123±85 

-1901209 

11128 

-99135 

Overall,  there  is  a tendency  to  overpredict  (iT  tends  to  be  nega- 
tive). This  is  az  might  be  expected  in  that  some  of  the  predictions 
were  made  conservatively  for  gage-ranging.  Also,  the  magnitude  of  n 
generally  increases  with  depth.  This  implies  that  the  magnitude  or 
the  overprediction,  at  least  relative  to  the  magnitude  of  the  measure- 
ment, increases  with  depth.  This  could  arise  because  of  failure  to  at- 
tenuate the  predictions  properly  (or  at  all)  with  depth;  or  it  could 
indicate  a difference  between  the  White  Sands  site  and  the  other  sites 
from  which  data  taken  served  as  the  basis  for  making  these  predictions 
by  scaling. 


D.  COMPARISONS  OF  THE  PREDICTIONS 


The  objective  of  this  section  is  to  compare  the  predictions  of 
AFWL,  WES  and  RDA  on  some  reasonable  basis.  The  sample  sizes  (sets) 
are  small,  sometimes  consisting  of  only  2 or  3 members,  so  it  is  not 
reasonable  to  make  critical  judgements  on  the  basis  of  individual  values 
of  a.  However,  the  question  of  how  the  three  techniques  fare  on  the 
average  may  be  asked.  We  now  form  a set  of  a for  each  agency  and 
again  compute  averages: 

N 

? = io,i 

where  only  positive  values  of  the  are  used  so  that  a must  always 
be  positive.  The  smaller  the  value  of  a,  the  better  the  consistency 
of  the  prediction;  that  is,  the  less  scattered  are  the  percent  devia- 
tion values  about  their  averages.  Compare  all  values  of  a from 
Table  7.3: 


TOTAL 


AFWL 

64 


WES 

84 


RDA 

60 


On  the  basis  of  this  comparison,  it  is  probably  not  possible  to  distin- 
guish between  the  AFWL  and  RDA  predictions,  but  the  WES  predictions 
rank  somewhat  below  the  others.  These  results  have  little  to  do  with 
conservative  overpredicting.  A negative  bias  in  n is  suppressed;  the 
comparisons  consider  only  predictive  scatter. 

Having  reduced  all  of  the  vertical  velocity  information  to  a single 
number  for  each  agency,  we  now  go  back  and  divide  the  numbers  in 
two  different  ways.  First,  we  consider  values  of  a calculated  separately 
for  the  direct-induced  (DI)  and  ai rblast-induced  (AI)  motions.  Second, 
we  compare  values  calculated  only  for  the  low  velocity  surface  layer 
with  those  calculated  for  the  underlying  higher  velocity  layer.  The 
results  appear  in  Table  7.4. 
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Table  7.4.  Values  of  a by  Motion  Source  and  by 
Seismic  Layer 


By  Wave  Motion  Source 


— 

AFWL 

WES 

RDA 

Direct- induced 

47 

84 

28 

Airblast-induced 

76 

- 

91 

By  Seismic  Layer 

Low  velocity  layer 
(0-3  m depth) 

30 

99 

51 

High  velocity 
substrate 
(>3  m depth) 

87 

72 

67 

For  these  cases,  the  picture  is  somewhat  unclear.  It  is  of  course 
possible  that  the  small  size  of  the  statistical  sample  may  be  affecting 
the  results  unduly.  Notice  that  prediction  consistencies  in  the  high 
velocity  substrate  show  a uniformity  among  the  agencies  which  is  not 
found  in  the  low  velocity  surface  layer. 

E.  CONCLUSIONS 

The  conclusions  reached  here  are  based  entirely  on  a study  of 
percent  deviations  of  the  vertical  components  of  particle  velocity. 

In  these  vertical  velocity  records,  the  airblast-induced  and  direct- 
induced  wave  components  can  be  distinguished  in  most  cases,  at  least 
in  the  slow  layer.  Tentative  conclusions  reached  in  this  study  may  be 
stated  as  follows: 

1.  All  three  of  the  agencies  made  reasonably  successful  predic- 
tions. This  is  probably  because  all  three  of  the  techniques  used  are 
based  essentially  on  scaling  from  previous  data.  Whether  or  not  the 
prediction  techniques  are  adequate,  and  how  they  should  be  improved, 
will  depend  on  the  needs  and  resources  of  the  weapons-effects  community. 

2.  Overall,  the  AFWL  and  RDA  techniques  seem  to  be  somewhat  more 
consistent  than  that  of  WES,  although  the  statistical  comparison  suffers 
from  inadequate  data. 
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The  conclusions  cannot  be  stated  more  strongly  because  there  are 
weaknesses  in  the  analytical  procedures  brought  about  by  the  data 
reduction  technique.  Additional  analysis,  beginning  from  a re-scanning 
of  the  original  data  traces,  will  be  necessary  if  a more  definitive 
result  is  to  be  obtained. 
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